Particle creation 
in the effective action method 

A.G. Mirzabekian^ and G.A. Vilkovisky^'^ 



^Lebedev Research Center in Physics, Leninsky Prospect 53, Moscow 117924, Russia 
^Lebedev Physics Institute, Academy of Sciences of Russia, Leninsky Prospect 53, Moscow 
117924, Russia 

Abstract 

The effect of particle creation by nonstationary external fields is considered as a 
radiation effect in the expectation-value spacetime. The energy of created massless 
particles is calculated as the vacuum contribution in the energy-momentum tensor 
of the expectation value of the metric. The calculation is carried out for an arbitrary 
quantum field coupled to all external fields entering the general second-order equa- 
tion. The result is obtained as a functional of the external fields. The paper gives a 
systematic derivation of this result on the basis of the nonlocal effective action. Al- 
though the derivation is quite involved and touches on many aspects of the theory, 
the result itself is remarkably simple. It brings the quantum problem of particle 
creation to the level of complexity of the classical radiation problem. For external 
fields like the electromagnetic or gravitational field there appears a quantity, the 
radiation moment, that governs both the classical radiation of waves and the quan- 
tum particle production. The vacuum radiation of an electrically charged source is 
considered as an example. The research is aimed at the problem of backreaction of 
the vacuum radiation. 
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1 Introduction 



Introduction proper. 

In the present paper we consider the problem of creation of particles from the vacuum 
by nonstationary external fields. This problem has been much discussed in the literature 
(see , e.g., [1-3]) but the present approach and the result obtained are new. Namely, we 
calculate the energy of massless particles created in external fields of arbitrary configu- 
ration and obtain the result as a functional of these fields. The calculation is done for a 
set of massless quantum fields coupled to all external fields entering the general second- 
order equation. The approach used is expectation-value theory and the method of the 
gravitational effective action. 

The main result of the present work is briefly reported in [4]. This final result is 
simple but its derivation is not because the effect of particle creation is nonlocal and sits 
in the cubic terms of the effective action. This is an effect of the one-loop vertices. A use 
of the nonlocal effective action in the context of expectation values is almost unknown. 
Therefore, we present here a systematic derivation showing the techniques involved. 

The result that obtains is remarkable since it brings the problem of the vacuum particle 
production to the level of complexity of the classical radiation problem. The strengths of 
external fields are expressed through their physical sources which are next integrated over 
certain spacelike hypersurfaces orthogonal to the geodesies. These integrals ( we call them 
radiation moments ) are direct generalizations of the moments of classical radiation theory. 
For external fields like the electromagnetic and gravitational fields, the radiation moments 
govern both the classical radiation of waves and the quantum particle production. As 
an example we consider the vacuum radiation produced by an electrically charged shell 
expanding in the self field. 

The work is motivated by the fact that the problem with external fields is physically 
incomplete. For a restoration of the energy conservation law it should be regarded as a 
part of a dynamical problem for the field's expectation values in an initial quantum state. 
An outstanding example is the gravitational collapse problem [5]. In the expectation- 
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value problem one first calculates the currents as functionals of the mean fields, and next 
makes these fields subject to the self-consistent equations. At the stage of the calculation 
of the currents the expectation-value fields appear as external fields. 

The setting of the problem with external fields in an asymptotically fiat spacetime 
assumes that these fields become stationary in the remote past and future in which case 
there exist the standard in- and out- vacuum states for the quantum field [1] . Because the 
external field is nonstationary in the intermediate domain, the out-vacuum is generally 
a many-particle in-state, i.e. the external field creates particles from the in-vacuum. 
The method usually applied for the calculation of this effect consists in representing the 
quantum field as a sum over modes and obtaining the Bogoliubov transformation that 
relates the basis functions of the in- and out- modes. 

The method of mode decomposition is normally used for an explicit solving of the 
equations in a given external field but one can also build a perturbation theory by solving 
for the basis functions iteratively in the external-field strength. This brings one to the 
loop diagrams of expectation- value theory, or the Schwinger-Keldysh diagrams [6-14]. As 
shown in [14], the Schwinger-Keldysh diagrams for the in-vacuum state are related through 
a certain set of rules to the Euclidean or Feynman effective action. This makes it possible 
to get away from both the mode decomposition and Schwinger-Keldysh diagrams, and 
carry out the calculation by the method of effective action. 

The vacuum energy-momentum tensor is obtained by varying the loops of the effective 
action with respect to the metric ( and next applying the rules of Ref. [14]). Therefore, 
irrespectively of the nature of the external field in question, one needs the effective action 
for the quantum field coupled to an external gravitational field ( in addition to the external 
field in question). Since, for the effect of the vacuum patricle production, the lowest 
nonvanishing order is second order in the field strength, one needs the terms in the effective 
action that arc qiiadratic in the field strength in question and linear in the gravitational 
field strength. The effect is thus contained in the one-loop triangular diagrams with at 
least one external gravitational line. 

The covariant perturbation theory for the effective action is built in [14-18] where 
the one-loop vertices are calculated for all couplings of the quantum field whose small 
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disturbances are propagated by the general second-order equation. These results make 
the starting point for the present work. 

Outline of the contents. 



In the present section, after a brief outline of the contents of the paper, we introduce 
the field model for which the calculation will be carried out, and the limitations under 
which the result will be obtained. We introduce also some notions pertaining to various 
limits at the asymptotically flat infinity, and relate the vacuum particle production to the 
radiation in the expectation- value spacetime. 

Sec. 2 reviews the structure of the effective action and the procedure of obtaining from 
it the expectation-value current which in the present case is the energy-momentum tensor 
for the mean metric. The energy-momentum tensor contains a contribution of the in- 
vacuum, , which is the subject of study in the paper. For the radiation problem, 
is needed only at the future null infinity and there are rules which facilitate 

obtaining T^^^ greatly. These rules are also reviewed including a theorem which relates 
the limit of X"*" for the kernel of an operator function to a certain limit for the function 
itself. This presentation is based on an earlier work [14-23] with the exception of the final 
result for 

The technique of this simplification is outlined in Appendix A. 

Sec. 2 concludes with a discussion of the problem of averaging the quantum noise 
which is a sign-indefinite contribution in T^^^ due to the quantum uncertainty. In two 
and four dimensions the mechanisms of this averaging are completely different since, in 
four dimensions, has the form of a total derivative in retarded time. The energy 

radiated by the vacuum for the whole history is thus determined by the limits of late 
time and early time, and, for this energy to be nonvanishing, certain nonlocal functions of 
external fields should exhibit a growth at late time. This offers a problem whose solution 
is the principal achievement of the present work. 

The result for T^*^ is obtained in Sec. 2 in the form of a superposition of nonlocal 
operators acting on the sources of external fields. All nonlocal operators in T^^^ are defined 



which is obtained by a drastic simplification of the results in Ref.[18]. 



by their spectral forms in which the resolvent is the retarded Green function [14]. A 
derivation of the kernels of these operators is the subject of Sec. 3. The resolvent is taken 
in the approximation in which it is determined by the geometrical two-point functions 
which, in their turn, are built on the basis of causal geodesies of the expectation- value (or 
background) spacetime. Therefore, all kernels in T^*^ inherit this geodetic structure and 
the retardation property. The geometrical objects involved are the light cone of a point, 
and the hyperboloid of equal (timelike) geodetic distance from a point. Subsequently, two 
other objects are derived from these: null and spacelike hyperplanes. 

The behaviours of all kernels at X+ are obtained in Sec.4. In particular, it is shown 
that, at the limit of X+, the kernel of the operator log(— □) which stands for the external 
line in T^^^ boils down to the kernel of !/□ . Here, there first appear the radiation 
moments but only their ultrarelativistic limiting cases defined as integrals over the null 
hyperplanes. Upon the calculation of the behaviours of the vertex functions at X+, there 
emerges a problem. It turns out that, in the case of the external vector field (and only 
in this case), the superposition of kernels in T^^^ fails to converge at T"*". The problem 
removes if the vector field is exposed to a special limitation that it contains no outgoing 
wave. The result in the paper is obtained under this limitation. 

As mentioned above, the total energy of vacuum radiation is determined only by the 
hmits of late time and early time. The purpose of Sec. 5 is a proof that the limit of early 
time makes no contribution. The proof uses three facts: i) the retardation of all kernels in 
, ii) the presence of a time derivative in the kernel of log(— □) , and iii) the assumed 
stationarity of the external fields in the past. Technically, it involves the properties of the 
geometric two-point functions and retarded kernels in presence of the Killing vector. The 
main result of this section is formulated as an assertion about the causality of the vacuum 
radiation. Sec. 5 discusses also the question of convergence of the masslcss operators l/D" 
for n > 1. This question emerges in connection with some of the vertex operators in T^^^. 

In Sees. 6 and 7 we come to the heart of the matter: the behaviours at late time. Most 
of Sec. 6 is devoted to the analysis of the behaviour of the function (l/U^X in the future 
of T^. The result is that the dominant contribution to this behaviour comes from X at 
the limit which here is defined as the limit of infinite proper time along the timelike 
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geodesies that reach the future asymptotically fiat infinity. It follows from this result that 
the radiation energy is given by an integral over the energies of the particles at i"*" and is 
nonvanishing only if the vertex functions have an appropriate growth at i'^. 

The behaviours of the vertex functions at i+ is the chief thing. They are obtained in 
Sec. 7. Here, there appear the full radiation moments defined as integrals over the spacelike 
hyperplanes. Much work in Sec. 7 is connected with the vector and tensor vertices. An 
alternative way of handling these vertices (which also requires much work) is considered 
in Appendix B. 

The final result for the energy of vacuum radiation is presented in Sec. 8, and it is 
gratifying to see its positivity. The positivity is based on two facts: i) the conservation of 
the currents of external fields, and ii) the sclf-adjointness of the equation of the quantum 
field. The discussion of the radiation moments is completed in Sec.8 by reviewing their 
role in classical radiation theory and establishing their relation to the textbook multipole 
moments. After that, it becomes visual that the quantum problem of particle creation is 
made almost the same thing as the classical problem of radiation of waves. 

In Sec. 9 the result for the vacuum radiation is specialized to spherical symmetry and to 
the external electromagnetic field. The longitudinal projections of the radiation moments 
are calculated. The vacuum radiation produced by an electrically charged spherical shell 
is considered as an example. The shell is assumed expanding in the self field from the 
state of maximum contraction to infinity, and the loss of its energy for the whole time 
of expansion is calculated in both the nonrelativistic and ultrarelativistic cases. Without 
accounting for the vacuum backreaction, the radiation of the ultrarelativistic shell violates 
the energy conservation law. The limits of validity of the technique based on the nonlocal 
expansion of the effective action are illustrated by specializing to the problem of particle 
creation in an external electric field. 

The concluding remarks are devoted to a discussion of the limitations imposed on the 
external fields. In particular, the emergence of the limitation on the vector field signals 
that the theory contains another effect: the vacuum screening or amplification of the 
electromagnetic waves emitted by a source. This effect is equivalent to an observable 
renormalization of all multipole moments and is analogous to the effect of the vacuum 
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gravitational waves [23]. 



The model. 

We shall consider the vacuum of a multicomponent quantum field = ({)^ coupled to 
external fields through the following equation for the small disturbances (J^?: 

m^p = Hj^dif^ = (1.1) 

with 

H = g^''V,VA + {P-^Ri) . (1.2) 

Here the hat over a symbol indicates that this symbol is a matrix in the space of field 
components: 1 = 6^^ , P = P\ , etc. The matrix trace will be denoted tr. In (1.2), P 
is an arbitrary matrix potential, R is the Ricci scalar of the metric g^^y , and is the 
covariant derivative with respect to any connection defining the commutator 

[v„ v.](5y. = = n^M^ . (1.3) 

The full set of external fields is thus a metric, a connection, and a potential, and the 
respective set of field strengths is 

3? = 7t^„ P) (1.4) 

where R^p^^, is the Riemann tensor of the metric in (1.2), and TZ^i, is the commutator 
curvature in (1.3). 

By solving the Jacobi and Bianchi identities 

v^7^„;3 + V/37^^, + v«7^^^ = , (1.5) 

Vcr-R^!//?" + Vy-Rcr^t/S" + ^fiRuais" = (1.6) 

with the vacuum initial conditions i, the commutator and Riemann curvatures can be 
expressed (nonlocally) through their contractions 

> = V,7^'^" , Jf"" = R^"" -^gf^'R (1.7) 

^We use the conventions for which [V^, S/^]X°' = R'pf.^X^^, R^p = ^1^/3, R = 9^'^Raf}- 
^i.e. with the condition of absence of the incoming waves. 
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(see Refs. [15-20,23] and Appendix A below). The set of quantities 

J = {Ji'\J>',P) (1.8) 

represents the physical sources of external fields (1.4). Specifically, the commutator cur- 
vature TZ^y is a generalization of the Maxwell tensor, and in (1.7) is a counterpart 

of the electromagnetic current. The expressions for the field strengths in terms of the 
currents are obtained iteratively but the conservation of the vector and tensor currents 
holds exactly: 

V^> = , . (1.9) 

Specifically, for the commutator curvature one has [24] 

[v^, v.]7^"^ = n^r.n'^^ - n'^^n^. + R'^^.n^^ - (i.io) 

whence 

v;\/,n^'' = ^ [v^, v.]7^'^'' = o . (i.ii) 

The vacuum energy of the quantum field is calculated below as a functional of the 
external- field strengths (1.4) to the lowest nonvanishing order in the number of loops and 
the power of 3?. Since solving the identities (1.5), (1.6) with the vacuum initial conditions 
is a part of the calculational procedure in covariant perturbation theory [15-18], the result 
will be expressed through the sources (1.8). At intermediate stages the potential term in 
(1.2) 

Q = P-^Ri (1.12) 

will as a whole figure in the capacity of a source but the final result will be expressed in 
terms of P for the reason explained in Sec. 8. 

The calculation in the paper is carried out under a number of limitations on the 
external fields whose significance is discussed in conclusion. One limitation is already 
predetermined: we consider only the fields of sources, i.e. neither the gravitational field 
nor the field represented by the commutator curvature will contain an incoming wave. It 
will be assumed that the sources of external fields have their supports in a spacetime tube 
with compact spatial sections and a timelike boundary. Their domain of nonstationarity 
will be assumed compact in both space and time. It will be assumed that the metric 
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has no singularities and horizons since here we consider only the case of finite particle 
production i. Finally, in the case of the commutator curvature there will be one futher 
limitation, namely that the vector source in (1.8) does not radiate classically. In the 
specific case where 7^^,^ is the Maxwell tensor, this limitation means that the external 
electromagnetic field contains no outgoing wave. No such limitation emerges in the case 
of the gravitational field. 

The limits X+ and i+ . 

An important role in what follows is played by the limits along the null and timelike 
geodesies traced towards the future. By the assumption above, all causal geodesies reach 
the future null infinity or the future timelike infinity (see, e.g., [25]). Here we introduce 
the notations and reference equations pertaining to these limits that will next be used 
throughout the paper. 

When dealing with the null geodesies, it is useful to build a Bondi-Sachs type [26,27] 
frame by choosing an arbitrary timelike geodesic (referred to as the central geodesic) and 
drawing the family of the future light cones with vertices on this geodesic (Fig.l). Let 

u(a;) = const. , {Vuf = (1.13) 

(with Vu past directed) be the equation of this family, 47rr^(a;) be the area of a 2 - 
dimensional section of a given cone in the induced metric, and (pix) be a set of two 
coordinates labelling the null generators of a given cone. The (f){x) takes values on a 
2-sphere and satisfies the orthogonality condition (V0, Vu) = 0. Then u is the retarded 
time, and r is the luminosity distance along the light rays that cross the central geodesic 
(the radial light rays). 

One property of the Bondi-Sachs frame used in the paper is the fact that, if the point 
X is in the causal past of the point x, then u{x) < u{x). Indeed, let o be the point at 
which the central geodesic crosses the past light cone of x, and x be any point belonging 

•^i.e. finite energy production. The number of created massless particles may be infinite. 
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to this cone or its interior. The past hght cone of x hes entirely outside the hght cone of 
o as illustrated in Fig. 2. We have u{x) — u{o), and u{x) < u{o) since all points for which 
u is bigger than u{o) are inside the future light cone of o. 

When calculating only in the asymptotically flat domain of infinite luminosity distance, 
the compact domain may be cut out. The choice of the central geodesic is then immaterial, 
the retarded time is measured by an observer at infinity, and labels the points of the 
celestial 2-sphere S to which the radial light rays come as r — > oo. 

To every radial ray there corresponds a 2-parameter bundle of parallel rays that come 
simultaneously to the same point of the celestial sphere and are, therefore, indistinguish- 
able at infinity. (They generate a null hypcrplane discussed in Sec. 4 .) Since the radial 
rays make a 3-parameter family, there is all together a 5-parametcr family of light rays in 
spacetime but only three parameters are distinguishable in the principal approximation 
at infinity: u, (p. Changing the central geodesic signifies going over to the parallel rays. 

We shall denote I'^[u,(l),r — > oo], or T'^fit, 0], or I'^, the limit of infinite luminosity 
distance r along the null geodesic that, when traced to the future, comes at the instant 
u of retarded time to the point of the celestial 2-sphcrc S. By using the arbitrariness 
u — > f{u), the retarded time will be normalized to coincide with the proper time of a 
distant observer at rest that registers the outgoing rays. This normalization condition 
reads 

(Vm, Vr) 

and the asymptotic form of the metric in the Bondi-Sachs frame is 



= -1 , (1.14) 



= -du'^ - 2dudr + r'^{de'^ + sin^ Odip'^) (1.15) 



where {9, ip) — (p & S. The integral over the 2-sphere S (normalized to have the area Att) 
will be denoted 

2lT TT 

J d^s{(t)) ^ J dip J de sine . (1.16) 





To build a finite vector basis at T+ define the scalar 



2r + u + 0{^) (1.17) 
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and the complex null vector tangent to the sphere S 

= r(V«^ + i sin ^V«9?) + 0(r°) . (1.18) 

x+ 

The resultant null tetrad at 

eaifJ.) ^VaU, VaV, rua, m*^ (a* = 1,2, 3, 4) (1.19) 

with II labelling the vectors of the tetrad, and m* complex conjugate to m satisfies the 
orthonormality relations 

(Vii)^ = (V^;)^ = = {Vu, m) = (V^;, m) = , 

(Vk, Vv) = -2 , (m, m*) = 2 (1.20) 
and can be used to expand the metric as follows: 



\v^uVj,v + Vf^vVvu) + ]:{m^ml + m*m,,) . (1-21) 



Consider now a timelike geodesic, and let s be the proper time along this geodesic. As 
s — > oo, the particle moving along the geodesic will go out of the domain of nonstationarity 

of external fields with the energy E (per unit rest mass). Only the geodesies with E > 1 
that reach the asymptotically fiat infinity are relevant to the present discussion. We shall 
replace E with the boost parameter 

, 0<7<1 (1.22) 

and denote i+[7,0, s — cxd], or i"'"[7,0], or the limit s — > cxd along the geodesic that 
comes to infinity with a given value of 7 to a given point (j) of the celestial 2-sphere. 

A distinction of this case from the case of null geodesies (apart from 7 = 1 in the latter 
case) is that the timelike geodesies differing by translations including the time translations 
are indistinguishable at infinity. There are again three parameters that register at infinity 
but the parameters are now 7 and 0, and there is a 3-parameter congruence of the geodesies 
that come to infinity with one and the same values of 7 and (j). (This congruence will be 
considered in Sec.7 .) All together, there is a 6-parameter family of timehke geodesies 
that come to infinity. 
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The Bondi-Sachs frame can be used also at with the asymptotic form of the metric 
in (1.15). For the asymptotically flat metric at ^"'"(7 > 0) one may use the vector basis 
(1.21) or the related basis 



(1.23) 



with 



t 



^r + u + const. , = - Vr , Vv^Vt + Vr . (1.24) 

The parameters of a point at or X+ can be replaced by a unit direction vector at 



infinity whose spatial components in the Minkowski frame at 



- dt + 5ikd-x!-d-xl 



;i.25) 



will be denoted rii = ni{(f)), i = 1,2,3. With the Euler parametrization (p — {9,(p) in 
(1.15), 

ni (0) = sin ^ sin (/? , 712(0) = sin ^ cos </? , n^lcp) — cos 9 , (1-26) 

and 

x^ = rn*(0) , n* = 5%fe , n'uk = I . (1.27) 
The basis vectors in (1.23) have in the Minkowski frame (1-25) the components 



(1.28) 



The limits and T+ are related. As 7 — > 1 the geodesic at i+ approaches the null 
geodesic that comes to at late time. Therefore, for an analytic function X, the sequence 
of limits i'^ and 7 — > 1 should coincide with the future of X+: 



2+[u,(p,r^oo] 



(1.29) 



under the substitutions 



s^- , 7 = 1-- + C>(— ) 
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(1.30) 
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Vacuum radiation in the expectation-value spacetime. 



In the framework of the problem of particle creation by an external field, it is impossible 
to answer the question where the energy of the created real particles comes from. Clearly, 
the vacuum particle production is only a mechanism of the energy transfer. The energy 
comes ultimately from the external field but this answer assumes that the external field 
should stop being external. 

In the self-consistent setting of the problem, the external fields become expectation 
values evolving from an initial quantum state. Specifically in the presently considered 
model, </7^ is the full set of quantum fields. Some of these fields have nonvanishing 
expectation values represented by the set of curvatures in (1.4). The expectation- value 
equations are obtained by the rules of Ref. [14] from the action 

1 /■ 

— J dxg^^'^R + 5'source + 'S'vac (1-31) 

where 5'sourcc is the (renormalized) classical action for TZ^^ and P, and S'vac is the sum of 
all vacuum loops. The equations for the expectation value of the metric can be written 
in the form of the Einstein equations 

R'''-lrR = 87r{TZr.^ + TZ) , (1.32) 

^source -^V2 Sg,. ^''''''^ 

in which there appears some extra source , the energy-momentum tensor of the 
in- vacuum. This source is a subject of the calculation below. Eq. (1-32) should be 
supplemented with similar equations for TZ^i, and P ( and, of course, there should be 
equations for the higher-order correlation functions, also derivable from (1.31), which will 
not be discussed here). 

The energy conservation law in the expectation- value spacetime is a consequence of Eq. 
(1.32) with the asymptotically flat boundary conditions. It is expressed in the existence 
of a conserved ADM mass Madm which equals the total energy of all sources and waves 
in the initial state, and a non-conserved Bondi mass M{u) defined by the metric at X+ 



16 



and accounting for radiation [25]. The difference M(— oo) — M{u) is the energy radiated 
from an isolated system by the instant u of retarded time, and 



M(-oo) = Madm . (1.34) 

The equation of the energy balance, or the Bondi-Sachs equation [26-27,23], is of the 
form 



dM{u) 



du 



;i.35) 



If d ^ 

— J d^S{(p) — Cgrav(M,0) 

X+ [u,(?!>,r— >oo] 

and is to be solved with the initial condition (1.34). The first term in (1.35) is the energy 
flux of the outgoing gravitational waves {dC^j..^^/ du is the complex gravitational news 
function [26,27] defined by the metric at X+), and the remaining terms are the energy 
fluxes through X+ of all sources of the mean gravitational field including the vacuum. The 
news function also contains a contribution induced by the vacuum [23]. 

We are presently interested only in the contribution of T^^^ to the radiation flux. 
This contribution will be nonvanishing only when the fields solving the expectation-value 
equations create real massless particles from the vacuum. If the other contributions are 
absent @ : 

dM{u) 



du 



J d'S{<P) {-/T!^:y (1.36) 

X+ [M,(/i,r— >oo] 

then the total radiated energy equals the total energy of created particles: 

M(— oo) — M(oo) = j = ^ ^^'^ kmft'^outl ill vac > (1-37) 

— oo P 

(see, e.g., [22]). Here a^^i^a^^^ are the creation and annihilation operators for the quanta 
of the field (p in the out-vacuum, and ep is the energy in the out-mode p. 

The inference from Eq. (1.37) is that the massless particles created from the vacuum 
are radiated through the future null infinity of the expectation-value spacetime. As seen 
from the initial condition (1.34), this vacuum radiation takes its energy from the ADM 



"^As is the case, for example, if the radiation of waves, both gravitational and matter, is banned by 
the symmetry of the initial state. 
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mass of the expectation value of the metric. Since the ADM mass is conserved, it can be 
calculated on a spacehke hypersurface taken in the remote past. The T^^^ is a retarded 
functional of the mean fields, and it vanishes in the remote past ( see Sec. 5 below). 
Therefore, the ADM mass remains unaffected by quantum corrections and equal to the 
energy in Tg'^uj.ce calculated on an initial hypersurface ( plus the energy of an incoming 
gravitational wave if any) i. It thus equals the energy of the "external" fields in the past 
from their domain of nonstationarity. In this way the energy conservation is restored. 

Eq. (1.36) with T^^^ obtained from the effective action will be used below to calculate 
the quantity (1.37). Special attention will be payed to the positivity of this quantity. 



^This has an important consequence that the ADM mass of the mean metric remains positive if T^q 
is energy-dominant [16]. 
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2 The vacuum energy- momentum tensor 
The effective action. 



The one-loop vacuum action for the field <^ is 

S^,,^'-Tr\ogH (2.1) 

where Tr denotes the functional trace, and H is the operator (1.2). Upon the calculation 
of the loop [14-18], this action takes the form 

-5va. = 5(2) + S(3) + 0[r] , (2.2) 
S{2) = dxg'/Hr ^^,{-D,)^,^,{i) , (2.3) 

□ = g^'-V^V, (2.5) 

which is the general form of a functional of the field strengths (1.4) expanded over a basis 
of nonlocal invariants [20]. The term S'(2) is a linear combination of five basis invariants 
5^i5^2(0 of second order in K with the operator coefficients (form factors) 7j(— □). The 
term S'(3) is a linear combination of twenty nine basis invariants 3?i3?23?3(i) of third order 
in 3? with the form factors ri(— Di, — Ds). A calculation of the effective action in 
any specific model or approximation boils down to obtaining the form factors in (2.2). 
The expansion (2.2) can also be assumed as a basis of a phenomenological theory of the 
vacuum [19-22]. 

The basis invariants of second order in 3fJ and their respective one-loop form factors 
are of the form [15] 

(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



3?i3?2(l) 


— p 




— RiR2^ , 




= PIR2 , 


3fJi3fi2(4) 


= P1P2 , 


3f?i3f?2(5) 
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7.(-n) - I^(>og(-^)-S) . (2.12) 

73(-°) = --^ , (2.13) 

74(-D) = -^log(-^) . (2.14) 

7.,(-n) = (2.15) 

where the parameter > accumulates the ultraviolet arbitrariness. This parameter 
doesn't affect the vacuum current at null infinity (see [21] and the calculation below). 

The third-order action contains no arbitrary parameters. Among the basis invariants 
of third order in 3fJ, eleven contain no derivatives, for example 

M23f?3(l) = AAA , (2.16) 
fourteen contain two derivatives, for example 

3?i3fJ23?3(12) = 7^l"^V^7^2;.aV'^7^3./3 , (2.17) 
three contain four derivatives, for example 

3fJi3?23?3(26) = V,V/3i?rV^V.i?2"^A , (2.18) 
and one contains six derivatives : 

5Ri5R2JJ3(29) = VxV^RfV^VpRr'^^V.Rti ■ (2.19) 

The full table of third-order invariants and their one- loop form factors is given in [17]. 

Since the operator arguments of the form factors in (2.4) commute, the form factors 
themselves are ordinary functions i . As analytic functions, they are defined by their 
spectral forms (in each argument) with the resolvents l/(n — m^). 

The action in the form (2.2)-(2.4) (i.e. with the loop done) determines both the matrix 
elements between the in- and out- vacua and the expectation values in the in- vacuum [14]. 
The difference is in the boundary conditions for the resolvents of the nonlocal operators. 



^Beyond third order in 3? this is no more the case [19,20]. 
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When the action (2.2) is varied, the resolvents are regarded as obeying the variational 
rule 

" ' ' ^^TT^. . (2-20) 



□ — w? □ — m? □ — m? 
and, after the variation has been completed, they are identified with the Feynman Green 

functions in the case of matrix elements and with the retarded Green functions in the 
case of expectation values [14]. Introducing a notation for the result of this procedure, 
we may write for T^^ in (1.32) the expression 

2 5S^ 



rpfXV 



^vac 



(2.21) 



ret 



For obtaining from the action (2.2) one needs to know the variational derivatives 
of the commutator curvature and potential with respect to the metric. The dependence 
of iZ^u and P on the metric is different in different models [24]. The calculation below 
is carried out for the case where Iz^y and P are metric independent but, as explained in 
conclusion, the final result is valid for arbitrary local TZ^i, and P . 

Obtaining T^^ at J+ . 



The T^^{x) is here to be calculated only for x at T+, and only the terms that con- 
tribute to the energy flux through are to be retained. These are the terms of order 



is used for 



X ^ with r the luminosity distance. Below, the notation T}^^ 
calculated up to terms 0(l/r^), and the terms 0{l/r^) are referred to as vanishing 
at 2:+. 

When computing TJ^^^ from the action (2.2) it is useful to have at hand the following 
behaviours that we quote from Refs. [21-23] and Sec. 4 below: 



□ 



= 0(1) , log(-n)xL^=o(l) , 



log(ni/n2) 



:X , logf-n)X =0(4) , (2.22) 

=0(4) . (2.23) 

The behaviours of these and other form factors including the higher-order ones can be 
obtained from the behaviour of the resolvent operator [21,22]: 



1 



□ 



X 



oc - exp(— Iconst.lmv^) 
x+ r 



(2.24) 
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The behaviours of the field strengths (1.4) and their sources (1.8) are in the most general 
case (see, e.g., Sec. 8 below) 



3? 



(2.25) 



and, with the present assumptions, these behaviours soften. Hence it is seen that the 
terms of T^^{x) in which the curvature K appears at the observation point x can often 
be discarded at X+. Also the terms in which there forms a positive integer power of the 
operator □ acting at the observation point can often be discarded since 



ax 



1+ 



Another useful fact is 



o{lx,x,) 



(2.26) 



(2.27) 



Both (2.26) and (2.27) follow from the form of the metric at X+, Eq. (1.15). 

Using Eq. (2.23) it has been shown in Ref. [23] that the variations of second-order form 
factors 7 make a vanishing contribution at X+. This fact makes an essential difference 
with the case of two dimensions (see below). Owing to this fact, the only contribution of 
the action ^"(2) to T^^ that survives at X"*" comes from varying the curvatures 3? in the 
products ^i^2{i)- 

Since varying destroys the curvature, the effective action of third order in 3fJ enables one 
to obtain T^*^ only to second order in 3?. Therefore, in the action S{3) only the curvatures 
3? in the products 3?i3?23?3(i) need to be varied. The contributions of variations of the 
third-order form factors F are already 0[3?^]. It is also useful to keep in mind that any 
commutator of the derivatives V with each other or with the form factors is 0[3?]. 

Using the rules above to discard the terms vanishing at X+ we obtain 



VcLC 



= ^trV^V^(7i(-n)+272(-n))m 

-^vnv^7l(-□)]i^-^v1v^7l(-□)]it:}i 



+ 
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_2_ 5S{?>) 



+ 0[3fJ^ 



(2.28) 
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where the variation of the action S{3) is not yet calculated but there appear the commu- 
tators which are of the same order as the variation of <S'(3). The commutators emerge 
when the terms linear in 3? are brought to their form in (2.28) with the aid of the Bianchi 
identities. Taking these commutators into account is necessary for to have the correct 
form at X+ (Eq. (2.36) below). Both commuting and varying of the operator functions is 
accomplished with the aid of their spectral forms [17,23]. By (2.21), all nonlocal operators 
in (2.28) have the retarded boundary conditions. 

Consider now the variation of the action S{3). A generic term in 5S{3) is of the form 



(2.31) 

X2=X3=X 



dxg'/^r{-0,, -n^, -n^)5U, ^2 ^3 + 0[3fJ3] (2.29) 
where has the structure 

6^{x) = VV6g{x) + Om . (2.30) 

Since, when going over to the variational derivative, the operator acting on dg{x) trans- 
poses, the respective term in the variational derivative at a point x has the following 
structure: 

vvr(-n, -02, -Ds)^{x2Mxs) 

where the operators VV and the first argument — □ of the form factor act at the obser- 
vation point X. In (2.31), first ^2 acts on ^{x2) and acts on ^{x^) with subsequently 
making the points X2 and x^ coincident with the observation point x, and next the argu- 
ment — □ of the form factor and the derivatives VV act on the thus obtained function of 
the observation point. 

The important fact is that one of the operator arguments of F, the one that in (2.29) 
acts on (53?, becomes an overall operator in (2.31). Therefore, the theorem in [21,22] 
applies by which the behaviour of the function (2.31) as x — > is determined by the 
behaviour of the function r(— □, —Da) as — □ — > with — □ the argument that 
becomes an overall operator in (2.31). For the function (2.31) to be 0(l/r^), the form 
factor should be [21] 

F(-n,-n2,-n3) = o(log(-n)) , -□ ^ . (2.32) 
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The log(— □) terms of the form factors at small □ determine the terms of at 

X+ [21]. Eq. (2.22) anticipates this fact. 

Thus, for the calculation of T>f^ at T+, one doesn't need the exact form factors. 
One needs only their asymptotic behaviours in each argument at the limit where this 
argument is small, and the remaining arguments are fixed. These asymptotic behaviours 
are presented in [18]. All the form factors F behave as follows {k = 1, 2, 3): 

r,(-ni,-n2,-n3) = i-AA(n„, □„) + log(-nfc)5>(n^, □„) + 0(nO) , 

—Of. —>■ , m ^ k, n ^ k, m < n (2.33) 

where Af' and Bj^ are functions of the two arguments □ other than Dfe. These behaviours 
differ from (2.32) since, in addition to the expected log(— □) terms, they contain senior 
1/n terms. However, as shown in [23], the l/D terms exactly cancel in the energy- 
momentum tensor. The variational derivative of S{3) is a sum of contributions of the 
form (2.31), and the form factors F enter this sum only in certain linear combinations. 
The 1/n terms cancel in these combinations leaving the log(— □) terms as the leading 
asymptotic terms at □ ^ [23]. 

Thus we infer that the contributions of the form factors F to are completely 

determined by the coefficient in (2.33). The table of these coefficients is given in [18]. 
The simplest ones are the coefficients in the form factor Fi of the basis invariant (2.16). 
They are of the form 

i5.'C„,D„) = iM5^E^ (2.34) 

for all the three values of k. The remaining B^^ are obtained by differentiating the 
generating expression 

1 □ — 1 □ \ ■ I 

Jm'—'m Jn'—'n 

with respect to the auxiliary variables jm,jn and subsequently setting j^n — jn — ^ 
[17,18]. All B/' are the thus obtained functions having also rational coefficients of the 
form or and rational additions. 

The table of the functions Bj^ in [18] can be considerably simplified with the aid of the 
technique outlined in Appendix A. The use of this technique makes it possible to bring 



vac 



to the final form below. 
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The result for T'^^ at X+ . 

vac 



The result for TJ^^ 



is of the form 

1+ 



1+ (47r)2 



^ V^Vlogi-n) I{x) (2.36) 



and is determined by a single scalar I{x): 

I{x) = trJi(x) + trJ2(x) + 0[SR^] (2.37) 

where ii{x) and /2(a;) are the contributions of first and second order in 3?. The contribution 
of first order obtains directly from (2.28) and from the expressions (2.11)-(2.12) for the 
second-order form factors: 

h = ~Ri . (2.38) 

The contribution of second order obtains from the results in [18] with the aid of the 
technique of Appendix A. It is of the following form: 

I2 = v^v^v^v, v^-^fr + v„v^v;:ft + v^v^K^L + Kcaw + trees , (2.39) 

= ^(^F(3,3)-lF(2,2))jrJ2"^i , (2.40) 
Kft = -^(Jf(2,2)-1f(1,1))jM" , (2.41) 
V:L = -^^{If{2,2)-^-F{1,1))j-^Q, , (2.42) 

V;calar = (^F(l, 1) - ^F(0, O))^^ (2.43) 

where J"'^, >, Q are the sources of external fields in (1.8) and (1.12), and F{n,n) with 
n = 0, 1, 2, 3 are specific cases of the vertex operator 



F{m,n) = {—) (— ) 



d \m/ d Nn log(jiai/j2D2) 



(2.44) 

ii=i2=i 



The V 's in (2.39) will be referred to as the vertex terms. As distinct from the V 's, the 
trees in (2.39) are the terms whose operator coefficients factorize into functions of one 



^The name "tree" is used here not quite in the usual sense. It should also be clear that both the 
vertices and trees are terms of the calculated loop. 
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variable: 

TREES= — fVfl-J"^)fV«-Jl)i + — J"^f-Ja/3)i-— ^f-i?)i • (2.45) 

The fact that only the sources of external fields appear in is a result of the use of 
the Jacobi and Bianchi identities in the calculation of the nonlocal effective action [15-20] . 
As will be seen, the conservation laws (1.9) for the vector and tensor sources play a crucial 
role in the consistency of the result above. 

The present calculation reveals a number of nontrivial properties of the effective action 
such as Eq. (2.36). When working with the exphcit one-loop form factors these properties 
emerge as a result of mysterious cancellations. In fact they can be predicted on the basis of 
the axiomatic approach to the effective action [19-22]. In particular, the behaviours (2.33) 
follow from the requirement of asymptotic flatness of the expectation- value spacetime and 
should hold for all form factors to all loop orders. The fact that the !/□ terms of these 
behaviours cancel in the energy-momentum tensor is predictable on the same grounds. 
As shown in [23], these terms stand for a vacuum generation of the gravitational waves. 
Eq. (2.36) is also a consequence of the asymptotic flatness. For the correct behaviour 
at X+ , T^^^ must have the form (2.36). Finally, the structure of the vertex terms in 
(2.39)-(2.44) should, as the results below suggest, be a consequence of the condition of 
unitarity encoded in the one-loop triangular diagrams. 

The problem of quantum noise. 



The vacuum energy-momentum tensor does not obey the dominant energy condi- 
tion. Therefore, the energy flux (^—dM/duj in (1.36) is not positive deflnite but the total 
radiation energy in (1.37) is positive. The point here is that the quantity in (1.36) is an 
expectation value over a quantum state rather than a classical observable. The indeflnite 
oscillations in (^—dM/duj should be within the quantum uncertainty, i.e. they represent 
a quantum noise which is present in even at the asymptotically flat inflnity X+ (see 
[22] for a detailed discussion). The problem is in discovering the mechanism by which the 
quantum noise gets averaged and the positive total energy emerges. 
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For the effective action in two dimensions 

5vac = const. J (fxg^/^R ^ R (2.46) 

this problem has a simple solution [5]. The positive radiation energy comes from the 
variation of the second-order form factor !/□ in (2.46). The contribution of this variation 
to is quadratic in R and energy-dominant whereas the contribution of the variation 
5R is linear in R and sign indefinite. However, this latter contribution has the form of a 
total derivative and vanishes in the full integral over time [5]. This is the mechanism by 
which the quantum noise sums to zero for the whole history. A counterpart of Eq. (2.46) 
in four dimensions is [15] 

Sy^^ ^ const. J d^xg^/^mog{-D)^ + 0[^^] , (2.47) 

and in this case the variation of the second-order form factor, log(— □), makes no con- 
tribution to T^^^ at X"*" [23]. The apparent problem in four dimensions is that is 
entirely a total derivative. Eq. (2.36). 

The insertion of (2.36) into (1.36) yields 

-'^=j^,fj£SW{r'lo,(-n)l) . (2.48) 

du [iny du^ J ^ ' x+[«,0,r^oo] 

Since this is a total derivative, the integrated energy flux is determined only by the limits 
of late time and early time u — ±oo. By the retardation property of the form factors , 
the contribution of early time vanishes (see Sec. 5 for the proof). There remains only the 
contribution of late time: 



M(-oo) - M(co) = Jim ^ A / cP5(^) (r^ log(-n)/ 



(2.49) 

X+ [ti,0,r— >oo] 



and, for it to be finite and nonvanishing, the integrand in (2.49) should have a linear 
growth at late time 



(rMog(-n)/) 



oc M , u ^ oo . (2.50) 

X+[u,(/),r— >oo] 



Since the setting of the problem with external fields assumes that these fields become 
asymptotically stationary as w — > ±oo [1], the growth in time required in (2.50) puzzles. 

27 



As we shall show, there can be no such growth if / in (2.50) is local in the field 
strength. The Ii in Eq. (2.38) is local, and, therefore, the vacuum radiation of first order 
in 3? is pure quantum noise [21,22]. The I2 in Eq. (2.39) is quadratic in 3? and nonlocal 
but the contribution of the trees is also pure quantum noise. The solution of the problem 
in four dimensions is that the growth in time is provided by the vertex operator (2.44). 
The kernel of this operator grows like u"'-+"--^^ ^ ^ the highest exponents m and 

n in all terms of (2.39) are precisely such that the result is (2.50). If the original quantum 
field contains no ghosts, the proportionahty coefficient in (2.50) is positive. 
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3 Retarded kernels of the nonlocal operators 
The retarded resolvent. 



With the form factors in the spectral forms, the only nonlocal operator in the current 
(2.21) is the resolvent l/[n — m^) which, in the case of the expectation- value equations, is 
the retarded Green function [14] . It admits an expansion of the same nature as the action 
(2.2), i.e. the expansion up to terms 0[3?"] by covariant perturbation theory. To lowest 
order in 3?, the retarded operator l/(n — m^) acting on an arbitrary tensor X^^ '^" is of 
the form 

past of X 

where a = a{x,x) is the world function [29,30] (one half of the square of the geodetic 
distance between x and x ), g^^ = g^fi{x,x) is the propagator of the geodetic parallel 
transport for a vector [30], Ji is the order- 1 Bessel function, and the integration point x 
is in the past of the observation point x. Here and below, the bar over a symbol means 
that this symbol refers to the point x. 

The geometrical two-point functions entering the Green function (3.1) satisfy the equa- 
tions [30] 

^''"V^dV.a = 2(7 , ^^^V^(7Vi,(7 = 2(7 , (3.2) 



a''V,g^^ = , a-^V,g"^ = , g% 



5\ , (3.3) 



x=x 



= -9% <^~' (3.4) 



with 



a^ = g^''a, , = , a, = , a^ = V^a . (3.5) 
To lowest order in 3? one can use the relations 

9% = -V^V-^a + 0[U] , (3.6) 

V.g% = 0[3?] , V^g%^0[lk] . (3.7) 
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As seen from (3.1), we are always dealing with some scalar source 

X = g%{x,x) . ..g%{x,x)X-^^--^-{x) (3.8) 

which may depend parametrically on the observation point. Since the resolvent will be 
used only in the approximation (3.1) i , we shall omit the symbol 0[3? x X] and use the 
short notation 

past of X 

which will always assume that the source on the right-hand side is transported in a parallel 
fashion to the observation point. 

For the integration over masses in the spectral integrals, the Green function should 
first be transformed by using the relation for the Bessel functions 



mJi{m^J—2(j) - - 



e^^MmV^) = ' (e-Va) (3.10) 

with an arbitrary timelike vector field C,^, and the Gauss theorem in the form valid for 
null boundaries: 

j dxg^'^V^A^^ = j dxg^/^A^^{y,T)5{j:) . (3.11) 

Here Vt is an integration domain, E = is the equation of its boundary, and S < holds 
inside Vt. Integrating by parts in (3.9) one obtains the expression 

-^^^X(x) = i- / dxf'H{-a)Urn^)V^{jJ^x) (3.12) 

past of X 

which holds with any choice of and from which the massless contribution proportional 
to 5((t) is absent. 

The kernels of !/□ , log(— □) , and of the vertex operators. 



The curvature correction to the lowest-order form factor is calculated in Sec. 4. 
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The retarded kernel of the massless operator !/□ is given by Eq. (3.9) with m = 0. 
Up to 0[3? X X] , 

-lx(x) = ^ J dxf/^6(a{x,x))x (3.13) 

past of X 

where the integration is over the past hght cone of the observation point x . 

The retarded kernel of the operator log(— □) has been calculated in [21,22]. To lowest 
order in 3? it is of the form 

1 r 

log(-n)X(a;) = — J dx^^/2(5'((T(a;,x))x+ a local term (3.14) 

past of X 

where the local term, i.e. a term proportional to X at the observation point x, is irrelevant 
to the present consideration i. 

The kernels of the vertex operators (2.44) are obtained by using the following spectral 
form of their generating function: 

log(jini/j2n2) _ 1 7 dm" 



ji^i - j2^2 3132 l {rn?/3i - / 32 - ^2) 

Combining (3.15) and (3.12) one finds 

log(jini/j2n2) ^ ^ , ^ 1 f -xi2^, I r 



XiX2(x) = ^ / (3.16) 

past of X 
past of X 

1 

— J dm'^Jo{m^J-2ai/ji)Jo{m^-2a2/j2) 



ji^i -32^2 ^ ^ 

past of X 



X 



where 

Gi = a{x, xi) , (72 = a{x, X2) . (3-17) 

The spectral-mass integral in (3.16) can be calculated with the aid of the Fourier-Bessel 
relation 

00 

/ dm^Jo{mJ-2ai/ji)Jo{mJ-2a2/32) = 26(— - — ) . (3.18) 

i 3i 32^ 



Expression (3.14) without a further specification is vahd only for x located outside the support of X . 

This includes the case where x X+ provided that X = 0{l/r^). The latter case is the one that we 

1+ 

presently consider. In the general case, the integral on the right-hand side of (3.14) is improper. For its 
precise definition and the form of the local term see [21,22]. 
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It is convenient to factorize the vertex function by introducing an auxiliary integration 
over a parameter. We write 

e[-a,)e{-a,)8{^ -^)- jdq 5{q - ^)s{q - ^) , (3.19) 

— oo 

and introduce the operator Hq depending on the parameter q and defined as follows: 

HqXix)^^ J dxf/'^S{a{x,x)-q)x , q<0 . (3.20) 

past of X 

Here the integration is over the past sheet of the hyperboloid a{x, x) — q associated with 
the point x (the past hyperboloid of x , Fig. 3). The operator Tig is a generahzation of 
the retarded operator — !/□ , Eq. (3.13). A calculation of the derivative 

d 



^^HgXix) = ^ J dxg^''^8'{q-a{x,x))X 



(3.21) 



past of X 



past of X 

with the result containing an arbitrary timelike vector field shows that this derivative 
with q replaced by jq will appear in (3.16) upon the use of (3.18) and (3.19). 
In this way we obtain 



d 



—— —— X^X2{x) ^2 dq ——Hgj.Xiix) 



d 



d{qj 



-Hq^^X^ix) . (3.22) 



Finally, the result for the kernel of the vertex operator (2.44) is 



^dq-' ^ '-^dq 



F{m,n)X,X2{x) = -2 I dqq"'+''[(-^y''^\x,{x)][(4:Y^\x2{x) 



(3.23) 



The kernels of the vertex operators superposed with !/□ . 



There are two more types of the vertex operators in (2.39): 

(a) ^^F{m,n) , (b) ^F(m,n) . (3.24) 
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It is important that in (a) n > 1, and in (6) m > 1, n > 1, i.e. the appearance 
of a factor l/Di or I/O2 in the vertex is necessarily accompanied by the appearance 
of a derivative d/dj bearing the same number 0. Therefore, it suffices to consider the 
generating expressions 



00 



1 a logfea./,,D,)^y^^^^ , ,.. Jf,.. „ ,, .Y..Y, . (3.25) 



1^2 di2 JiDi - 32^2 jij2^ / {.rn?/3i - ^i){jn^ / 32 - ^2)^ 







00 



1 d d log(jini/j2n2) ^ v - ^ [ ^"^^ v v 



0^^2831832 31^1-32^2 31^32^ { (mVji-ni)2(mVj2-n2)2" 

(3.26) 

The square (and, generally, any power) of the massive Green function is a well defined 
operator whose kernel can be obtained by differentiating the Green function with respect 
to the mass: 

1 8 1 

>X = -^,j-^—,X . (3.27) 



(m^ — n)^ 8m? {rn? — □) 

The behaviour of the function (3.27) as m? and hence the convergence of the integrals 
(3. 25), (3. 26) at the lower limits depends on the class of the sources X. This question will 
be considered in Sec. 5. 

The differentiation of the Green function with respect to the mass is readily accom- 
plished in the original expression (3.9) by using the relation for the Bessel functions 

8 mJi(mV —2a) 1 



- Jo{m^/^) . (3.28) 



8m^ v/^2^ 2 
The result is 

-^^—-X{x) = ^ I dxf/'0{-cx)Jo{mV^)X . (3.29) 

past of X 

The insertion of (3.29) and (3.12) into (3.25), (3.26) leads to the same spectral in- 
tegral and the same factorization procedure as in the previous case. For the generating 
expressions we obtain 

1 8 log(jini/j2n2) Y vf \ } ^ \ ^ %J V I \\\'^ %J V i \^ QnA 
-— T^- — ——X^X2[x)= dq ——Hgj.Xiix) -ngj,X2{x) , 3.30 



^"This is a manifestation of the general property of the form factors estabhshed in [17] with the aid of 
the integral a-representation and called there " rule of the like a" . 
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_l_ d d log(jiDi/j2a2) 



XiX2{x) = - dq —Hqj^Xi{x) —Hqj^X2{x) 



32 



(3.31) 



and, finally, the results for the kernels of the vertex operators are 



^F{m,n)X^X2{x) 

^2 



(3.32) 



d \n~l 1 



-/ '-n,x,{x) , n>l. 



F{m,n)X^X2{x) 



d \m-l 1 



d \n-l 1 



dq' 



X] 



(3.33) 
, m > 1, n > 1. 



The question of convergence of the integrals (3. 25), (3. 26) at = transfers now to the 
integrals (3. 32), (3. 33) with respect to q whose convergence at ^ = — oo should be provided 
by the properties of the sources X. 
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4 The asymptotic behaviours at X+ 



Null hyperplanes. 



The power of growth of the world function as one of its points tends to and the 
other one stays in a compact domain should be the same as in flat spacetime. Therefore, 

we may write 

a{x, x) =rZ{u,(j); x) + 0{r°) (4.1) 

a;6X+[u,</>,i — >oo] 

with some coefficient function Z. 

Inserting the behaviour (4.1) into the equation (3.2) for a with respect to the point x 
and using the asymptotic form of the metric at X+ we obtain 

d 

-2rZ—Z + 0(r^)=2rZ + 0(r^) (4.2) 
ou 

whence 

^Z{u,<P; x) = -l . (4.3) 
Therefore, we introduce a new notation to write down the solution for Z 

Z(u, ; x) = -li + U^(x) , (4.4) 

and rewrite Eq. (4.1) in its final form 

a{x, x) ^r(-u + U^{x)) +0{r^) . (4.5) 

Inserting the behaviour (4.5) into the equation (3.2) for a with respect to the point x 
we obtain 

(W(x))' = , U{x) = U^{x) . (4.6) 
Hence we infer that the equation 

U^{x) = u (4.7) 

with fixed cj) defines a family of null hypersurfaces labelled by the retarded time u but 
different from the future light cones of the Bondi-Sachs frame. The hypersurfaces (4.7) will 



35 



be called hyperplanes . To different directions at infinity there correspond different 
families of null hyperplanes. 

An important property of the function U^{x) follows from Eq. (3.4). Inserting the 
behaviour (4.5) in (3.4) one obtains 

Vf,U^{x) = ^7/(x, x)V^u - (^/(x, x)V^r^U^{x) (4.8) 



where a = 1, 2, and x . By (1.1^ 



oc -m^ 

r 1+ 



(4.9) 



while the contraction gfi^{x^ x)m^ remains finite as a; X+ [23]. Therefore, 

g-^'^{x,x)V,r =0(-) . (4.10) 

As a result, one obtains the following law of parallel transport: 

V-^U^{x)=gJ'{x,x)V^u{x) . (4.11) 

The geometric meaning of the equations above will clarify if one answers the following 
question: what becomes of the past light cone of a point x at the limit x X+? At 
this limit, one of the null generators of the cone merges with i.e. disappears from any 
compact domain. Therefore, the resultant limiting surface is no more a cone although 
it remains a null hypersurface. As follows from Eq. (4.5), this limiting surface is none 
other than the hyperplane (4.7) whose parameters m, label the point of X+ to which 
the vertex of the prelimiting cone comes (Fig.4). The null generators of the prelimiting 
cone all but one become the generators of the limiting hyperplane. It follows that the 
generators of each hyperplane in (4.7) are the null geodesies that, when traced towards 
the future, come to one and the same point (f) of the celestial sphere at one and the same 
instant u of retarded time. This can be taken for a definition of light rays parallel in the 
future. Since the null generators of a hyperplane merge at infinity, they have a common 
tangent vector at 1^ which is Vm. On the other hand, at a point x of a compact domain 



^^More specifically, they should have been called future hyperplanes as distinct from the past hyper- 
planes defined similarly by the conditions at . However, the past hyperplanes do not figure in the 
present consideration. 
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the vector tangent to the generator is VpU^{x) . Eq. (4.11) is, therefore, the law of 
parallel transport of the tangent vector along the null geodesic emanating from a given 
point and belonging to a given hyperplane. 

Further properties of the function U^{x) derive under the assumption that for both 
the point x and the point at X+ one can use one and the same global Bondi-Sachs frame: 

U^{x) = U^{u,^,f) . (4.12) 

Then we have 

U4u,^,f)>u , (4.13) 

and 



U^{u,(f),r) 



d _ - _ 
u , -yf/0(u, 0, r) 



. (4.14) 



The inequality (4.13) is a consequence of the general fact mentioned in Sec.l that, for the 
points X belonging to the past light cone of x , u{x) < u{x) . With x at X+ , the latter 
inequality holds for the points x belonging to the limiting hyperplane (4.7). Eqs. (4.14) 
follow from the fact that the radial geodesic u = u , (p = (p along which the point x in 
(4.5) tends to X"*" serves at the same time as a generator of the past light cone of x and, 
therefore, belongs to the limiting hyperplane (4.7). The second equation in (4.14) implies 

Vt/(i _ = Vu and follows from the first equation and Eq. (4.6). 

(f>=(l> 

Finally, in the case where x is at the future asymptotically flat infinity X+ or one 
can use the flat-spacetime formula for U^{x). Indeed, in this case the geodesic connecting 
the points x and x in (4.5) passes entirely through the asymptotically fiat domain.0 One 
has 



UAx] 



= u + r{ i — cos uj 



+ f(l 

t-ni(0)x^ (4.15) 



where the first form refers to the Bondi-Sachs parametrization (4.12) and the second to 
the Minkowski coordinates (1.25) for x. In (4.15), uj{(j), (p) is the arc length between the 



^^The geometric two-point functions cr(x, x) and g^^{x, x) are nonlocal objects involving the metric 
on the geodesic connecting the two points. Therefore, even if both points are at the asymptotically flat 
infinity but one is in the future and the other in the past, the flat-spacetime approximation is invalid 
since the geodesic connecting the two points passes through a domain of strong field. 
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points (j) and on the unit 2-sphere: 

cosuj{(f), (j)) = ni{(f))n\(f)) 
Here ni{(p) is the direction vector introduced in (1.26). 



(4.16) 



The !/□ and log(-n) at 1+ . 



The behaviours of the kernels of l/D and log(— □) at T+ follow immediately from 
(3.13), (3.14) and (4.5): 



□ 



= --Di{u,<l>\X) , (4.17) 

X+ [«,</), r—>c»] 1" 

= -^|-Di(«,0|X) , (4.18) 



\og{-n)x 

and the coefficient of these behaviours is an integral over the null hyperplane (4.7): 

Di{u, (j)\X) = ^J dxg^'^5{U^{x) - u)X . (4.19) 

As discussed in Sec. 7 below, this integral is the ultrarelativistic limiting case of a more 
general hypersurfacc integral D called radiation moment of the source X. The subscript 
1 that D\ in (4.19) bears is to distinguish the ultrarelativistic moment. 

It follows from Eqs. (4.18) and (2.48) that the moment Di determines the flux of 
vacuum energy through T+. We obtain 

dM{u) 2 



j (fS{ct>)Dr{uA\I) (4.20) 



du (47r)2 dvP' 

where / is the scalar in (2.36). Since D\ serves at the same time as a coefficient in (4.17), 
it governs also the classical radiation (see, e.g.. Sec. 8 below). Of special importance 
is, therefore, the question of convergence of the integral (4.19). It is only with this 
convergence that the behaviours (4.17) and (4.18) hold. 

Note that the convergence of the original integral (3.13) defining the retarded Green 
function imposes a restriction only on the behaviour of the source X at the past null 
infinity whereas the convergence of the moment (4.19) requires certain behaviours of X 
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also at the future null infinity and spatial infinity. The behaviour at the future null infinity 
is critical. 

To study the convergence of the integral (4.19) at consider a portion of the inte- 
gration domain defined by the following inequalities in the Bondi-Sachs frame: 

Q : U2> u> Ul , 4> + A> 4>> 4>- A , f > ro (4-21) 

with large ro, small A, and any fixed ui, U2. In the integral restricted to the integration 
point X is near X+. Therefore, one can use the fiat-spacetime expressions for the metric 
and the function U^{x), Eqs. (1.15) and (4.15). One finds 

U2 (p+A 



J dxf/^d{u^{x) -u)x = JduJ d^S{(p) J dff^d{u-u + f{l-cosLu))x 

= /,/rj^(.-A-)^J.,^,f=-ii^] (4.22) 

J J 1 — COS LU^ ^2:+ 1 — COS UJ 

Ul </)— A 

with cu — uj{(t)i 4>) in (4.16), and f large i.e. (1 — cos a;) small. In agreement with the fact 
that the null hyperplane reaches the future infinity at only one value of , = , the 

question of convergence concerns the angle integral. The source X may have any power 
of decrease or even increase at X+ 



X 

provided that 



= k(u, (fyf^ (4.23) 

X+[u,</),f— >oo] 



= O((0- 0)5+2") . (4.24) 

However, to have 5 + 2n 7^ in (4.24), the source X should depend on the external angles 
(j). This is possible in the case of a tensor X since in this case X depends parametrically 
on the observation point through the propagators of parallel transport, Eq. (3.8). For 
a scalar source that has no special relation to the integration hyperplane in (4.19), the 
convergence condition is 5 + 2n = , or, with the analyticity of X at T+ , 

X =0(4) . (4.25) 
x+ ^ r'^ ^ 

The / in Eqs. (2.36)-(2.39) and (4.20) is a scalar. We must, therefore, make sure that 
it satisfies the criterion (4.25): 
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Let us check if this criterion is fulfilled for the vertex contributions. 



The vertex operators at . 



The asymptotic behaviours of the vertex functions (3.23) as x — > X+ are obtained by 
making the following replacement of the integration variable in (3.23): 



q — r{T — u) 



where r is the new integration variable, and r are the parameters of 



x e 0, r — > oo] 



With q replaced as in (4.27), and r fixed. 



^-D^{T,cj>\X) 

a;€X+[u,0,) — >oo] f 



As a result we obtain 

F(m,n)XiX2 
2 



X+[u,(/),r— >oo 

d \ m+1 



= / <;r(r-«)-»[(|:)'""c.(r.^|X0][(-i)-'D.(r,^|A-.) 



In the same way, for the vertex functions (3.32) and (3.33) we obtain 
^■F(m,n)XiX2 

X+\u,<l>,i — voo] 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



i / „r«[(A)"'«c,(r,^|A-0][(^)"-'-l-Z,,(r,^|A-.) 



-F(m,n)XiX2 



1^2 



2+[u,<j),i — >oo] 



(4.31) 



(r-u)' 



Thus we infer that, with the general sources, the operator in (4.29) docs have the 
needed power of decrease at X+ but the operators in (4.30) and (4.31) do not. In Eqs. 
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(2.39)-(2.42) these latter operators act only on the vector and tensor sources and are 
accompanied by extra derivatives but the derivatives alone do not help since they have 
the projections tangential to What helps is the conservation of the vector and tensor 
sources. Owing to this conservation, the tangential derivatives vanish, and the behaviours 
in (4.30) and (4.31) soften down by at least one power of 1/r. As will be seen below, in 
the case of the vector vertex this is still insufficient and, to ensure condition (4.26), we 
shall have to impose a special limitation on the vector source. 

The conserved charges. 



Below, for the moments Di of the sources in (1.8) we use the short notations 



Di{rA\P) 



D, (t. 



Di 



D 



flU 



(4.32) 
(4.33) 
(4.34) 



To distinguish the moments of the Ricci scalar R and the matrix Q in (1.12), they will 
be denoted 



Di(r,</.|g) = Di«(r,0)=Di« 



(4.35) 
(4.36) 



We have 

9,.Dr = -D^^ . (4.37) 
The non-scalar moments contain the propagators of parallel transport. Specifically 0, 

(4.38) 



x^X+ [r,( 



= dxflH[U^{x)-T)g%{x,x)P{x) 
Dr = ^/t/x^'/'5(f/<^(x)-r)(?^^(x,x)^?'^,(x,x)J^''(x) 



(4.39) 



X^I+lT,(f)\ 

and it is important that the propagators g°'^{x, x) are taken at the point x of X+ with the 
same parameters (t, 0) as in the 5— function. The moments D\^D-^,D-^^ , etc. have the 



^■^The matrix moments contain implicitly the propagators of parallel transport for the matrix indices. 
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transformation properties suggested by their indices but it should be remembered that 
they are defined only for a point at X+. Thus, the metric Q/j^i, in Eq. (4.37) can only be at 
the point (r, 0) of T+. Eq. (4.37) itself is a consequence of the relation [30] 

9,.u{x)g%{x,x)g''^{x,x) = gpt^ix) . (4.40) 

Using the law of parallel transport in (4.11) we may calculate 

VaU^Di'' = -^jdxf'''5'(y^{x)-T)Vj,U^{x)J~^{x) 

= -^j dxg^'^V-^5(u^{x)-T)P{x) . (4.41) 

Since, by assumption, the support of is confined to a spacetime tube, the integration 
by parts in (4.41) gives rise to no boundary term. Hence, using the conservation law (1.9) 
we obtain 

V^u^b^^^Q . (4.42) 
The conserved quantity in (4.42) 

e = V«'u£>i" (4.43) 

is the full charge of the vector source since the charge can be written as the integral 

1 f 

e = —j dxf^^S(E{x))VpE P{x) (4.44) 

over any hypersurface E = crossing the support tube of (VS is past directed). In 
expression (4.41) this hypersurface is a null hyperplane. If is a pure vector like in 
the case of the electromagnetic current, Eqs. (4.42)-(4.44) are exact. In the general case 
where the integral (4.44) contains the matrix propagators of parallel transport, they are 
valid with accuracy 0[3?^] since the proof of conservation makes use of Eq. (3.7). 

The significance of the assumption about the support of is in the fact that it 
excludes a radiation of the charge. In the general case, the full conserved charge is given 
by Eq. (4.44) in which the hypersurface E = is spacelike. With this hypersurface null, 
the charge is generally not conserved. It particular, if this hypersurface is the radial light 
cone u — const., Eq. (4.44) defines the Bondi charge 

1 f 

e{u) = ^ y dxg^/^5{u - u)Vp,u J^{x) . (4.45) 
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Then we have by using Eq. (3.11) : 

= - ^lim^ ^ j dxf/^S{u -ujS{v- vo)Vp,vJf'{x) (4.46) 

with V in (1.17). The boundary v — vq with — > oo is X+. Since at this boundary one 
can use the asymptotic form of the Bondi-Sachs metric (1.15), the result is 



X+fM,<i,r— »ool 



(4.47) 



which is an equation analogous to (1.35). If, on the other hand, the support of is 



confined to a tube, we have 



T+ 



, and de/du — . 



In the case of the tensor source, a calculation similar to (4.41)- (4.42) and vahd with 
accuracy since Eq. (3.7) has to be used yields 

d 



The respective conserved charge 



(4.48) 



(4.49) 



is the full momentum 



p'^ = i-/ dxg''\T.{x))g%{x,x)V-;>: ^'{x) 



(4.50) 



and a radiation of the charge is again excluded by the assumption about the support of 



The vector and cross vertices at . 



In Eq. (4.29) the moments Di always appear differentiated with respect to time. This 
is no more the case in Eqs. (4.30) and (4.31) but in the specific combination of F{m, n) 
entering the vertex functions (2.41) and (2.42) the terms with the undifferentiated D\ 
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cancel. Using integration by parts we obtain 



X+[u,(t>,r^oo] 



1 " 



d 



^^i"^(^''^)J[^^i''(^''^) (4-51) 



dT 



yap 



vect 



1 

12^ 



(4.52) 



X+ [u,(/),r— >oo] 

Owing to the presence of the derivatives acting on the Di's in all terms of these expres- 
sions, we have 



V a" I' cross 



x+ 



= , Vo^uV, 



al3 
vect 







(4.53) 



by (4.42) and (4.48). 

For a tensor V°'^ with the property (4.53) we are to calculate the behaviour at 
0,r — > oo] of the quantity Va^/sV"'^ appearing in (2.39). For that, we expand 
yaP Qygp ^jjg null-tetrad basis in (1.19): 

= e"(Af)e^(z/)l^(H (4.54) 

and use that 







by (4.53), and 



Ve"(/x) 



x+ 

by the result quoted in [23]. 

Since V{ijlu) is a collection of scalars, one has 



0(1) , VVe"(/.)L^=0(l) , etc. 



(4.55) 



(4.56) 



whence, in view of (4.9), 

Va^(H 

Similarly, 

VaV^V(H 



(4.57) 



(4.58) 



q2 q 



(4.59) 
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Upon the insertion of these expressions in VqV pV^^ and the use of (4.54)-(4.56) along 
with the condition V{e°'{^)Vau) — , one is left with 



X+ 



(4.60) 



An explicit calculation in the metric (1-15) yields 



1 



1+ 



whence finally 



By using (1.21) and (4.53), this result can also be written in the form 

1 d 



r du 



(4.61) 



(4.62) 



(4.63) 



One power of 1/r is thus gained. 

For y^fss (4-51) this gain is sufficient to satisfy the criterion (4.26): 



X+[i(,</),r— >oo] 



d 



(4.64) 



but in the case of V^ect i^i (4-52) we wind up with the behaviour 1/r^: 

(4.65) 

The integrand in the latter expression 



dr 



dr 



(4.66) 



is none other than the energy flux of the outgoing waves of the vector connection field 
(see Sec. 8). To ensure the fulfilment of the criterion (4.26) we are compelled to impose 
a limitation on the vector source, namely that this source does not radiate classically; 
then the quantity (4.66) vanishes. The energy of the vacuum radiation is obtained in the 
present paper under this limitation, and the limitation itself is discussed in conclusion. 
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The gravitational vertex at X+ . 



Using expression (4.31) for the F{m,n) entering the vertex function (2.40), and inte- 
grating by parts, we obtain 

= ~ / cLt^ [t — u}" - 

1+ [u,()>,r^oo\ 



grav 



360r 



(4.67) 



+(r-ur[^D^'^^{r,<P)] [-^Dnr,cl>)]+2D^'^^{r,<t>)Drir,<P)] ■ 

Here, as distinct from the previous case, the terms with the undifferentiated Di do not 
cancel. For the completely symmetrized V^i-^'^ we have 



1+ [u,(f>,r—fOo] 



^-p^^ ] dTDr\rA) (4.68) 



by (4.48) and (4.49). Two more time derivatives and two more contractions with Vu are 
needed for this quantity to vanish: 



rs2 



x+ 270r du 



V^U V/JU Vq;U ^gra.v 







(4.69) 
(4.70) 



Above, the tensor V^rav"' considered symmetrized because, with the appropriate 
accuracy, it is symmetrized in the quantity appearing in (2.39): 



(4.71) 



For the calculation of this quantity at I^, we introduce the projection of V^^£'^''^ on the 
null tetrad 

V^t""^ = e^^)e^{a)e'^{p)e^{6) V{^apS) (4.72) 

and use that by (4.70) 



rj2 

V^uVpuVaU e"(7)e^(a)e'^(p)— \/(7ap5) 







1+ 



(4.73) 
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One needs only the third and fourth derivatives of Vl'japS) since by (4.56) 
V«V^V^V.yg(:^'^'^)|^^ = 4v4e"(7)e'3(a)e^(p)e^(5))V;3V^V.V^(7ap<5) (4.74) 

The third derivative is needed only to lowest order, 

V{'yap6) + 0(^-V) 



(4.75) 



whereas of the fourth derivative one needs two orders: 



VaVpV^V^V{-fapS) 



= VauVf3uVi^uVuU-^V{-fap5) 

+ {y qUV I3UV ^uV yT + V aUV i}UV yUV 

+VauVfj,uVyuVpr + V puV ^juV yuV oct) q^zq^Y ^'^^P^^ 

+ {y aUV puV u,uV yCj)" + V aUV I3UV yUV ^(j)'' 

^-VyuVfjU ■ Vy/aU + VauVpu ■ V yUj —V{-fapS) 
+ 0{^,V) . (4.76) 



All terms in the two latter expressions have a sufficient number of d/du for Eq. (4.73) to 
work but not all have a sufficient number of Vu. As a result, one is left with 



' grav 



= -6V„«V^«(V,V.«)— V;(;:r^ + 0{-V) . (4.77) 

It is already seen that the present case is dissimilar to the previous one (cf. Eq. (4.60)). 

Since Wu = 0{l/r) , there can be and is only one factor of VVm in (4.77). The 
remaining indices have to be contracted with Vu. Therefore, the completely transverse 
projection of V'^f^'' drops out of (4.77). 
Using (4.61) we have 

X. ^ r rn.mlV^ir^) + 0{-V) , (4.78) 
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or, by (1.21) and (4.70), 

VoSpV.V.V-tr ^= -a^(Va«V^« 9,.V^fr^) + 0{-^v) . (4.79) 
Finally, using (4. 69), (4. 37), and denoting 

^ = V.iip" , (4.80) 

we obtain 

V.V,V,V.V',trL^^_^_=;^i^;^ *) + 0(i) . (4.81) 

The result is that the 1/r^ term survives also here but there is nothing like a radiation 
flux in its coefficient. The most important fact is that this coefficient vanishes at late 
time as distinct from the coefficient in (4.65) which grows at late time (Eqs. (6.53) and 
(6.54) below). Since only the late-time behaviour is relevant to the calculation of the 
total vacuum energy, the term in (4.81) will have no effect on this calculation and 
no special limitation on the tensor source will be required. 

The convergence condition (4.26) should, of course, be fulfilled everywhere at T+. The 
full solution of this problem including a removal of the limitation on the vector source 
remains beyond the scope of the present work (see the conclusion) . 

The trees at X+ . 



There remain to be considered the terms in (2.38)-(2.39) whose operator coefficients 
factorize into a product of massless Green functions. The simplest such terms are Ii in 
(2.38) and the trees in 12 that contain a curvature at the observation point : 

R , J'^^i^Je.fs) , R{^R) . (4.82) 
Generally, in an asymptotically flat spacetime one has (see, e.g., [23]) 




(4.83) 
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and the power of decrease of (l/n)X at X+ is 1/r. Therefore, in the general case all the 
trees in (4.82) are 0(l/r^), and with the present assumption about the support of the 
sources they vanish identically outside a spacetime tube. 
The remaining tree in (2.45) can be written in the form 

1 



□ 



where 



and use is made of the conservation law (1.9). By (4.17) 



and, to lowest order in 1/r, 



1+ 



VauVpu^A''^ + 0{-A) 



(AM) 



(4.85) 



(4.86) 



Hence 



1 

du 



^g^.p'-p" + O 



(4.87) 



(4.88) 



where p'^ is the conserved momentum (4.49). Thus, also this tree is 0(l/r^) at X+. 

There is, however, one more tree which thus far has not been taken into account. 
The point is that we are using everywhere the lowest-order approximation (3.1) for the 
resolvent. This is correct in the terms of second order in 3? but in the terms of first order 
in 3? the resolvent itself should be taken with a higher accuracy. Therefore, we must come 
back to expression (2.36) with the first-order J , / = Ji , and calculate the curvature 
correction to the form factor in this expression: 



x+ 



(51og(-n)Ji 

By using the spectral form 

oo 

log{-n/c') = -J 



^*The curvature corrections to the derivatives V^V^ in (2.36) vanish at X+. 



(4.89) 



dm' 



1 



2 — □ 77^2 _|_ ^2 



(4.90) 
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and the variational law for the retarded Green function we obtain 

51og(-n)X(a;) = -Jdm^^^^^m,^^^^X{x) = ^^^^^6n,X{x) (4.91) 

where the numbers on the D's indicate the order in which the operators act on X{x). We 
may now use the theorem in [21,22] by which the limit of (4.91) as a; —> X"*" is determined 
by the limit □ — > in the operator □ acting the last i.e. in O^. We obtain 

5log(-n)X ^ - ^^^^~°^h n2X ^ -log(-a)(5D-X . (4.92) 

As could be expected, the variation of log(— □) in the background fields behaves at X+ in 
the same way as log(— □) itself, i.e. like 1/r^, and the result (4.92) is valid up to 0{l/r^). 
Using the latter result in expression (4.89), we may write 

(log(-n) + 51og(-n))/i|^^= log(-n)(7i - 5n^i,) . (4.93) 

This means that introducing a correction to the form factor boils down to the replacement 
of the scalar I hy I + SI with 

51 = -50^Ii . (4.94) 

Upon this replacement one may use everywhere the lowest-order approximation for the 
resolvent. In (4.94) one is to insert the expression for h from Eq. (2.38) and the expression 

SO = 2(ii?"^)VaV;3 (4.95) 

from Ref. [15]. As a result, the correction under the sign of trace in (2.37) takes the form 

or, in terms of the conserved current J"^, 

« = i;V„V,B"^ + ^fl(ifl) (4.97) 

with 

^ {^r^){^R) . (4.98) 

Expression (4.97) is one more tree. 
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The second term in (4.97) is of the type (4.82), and its contribution at X"*" is 0(l/r^). 
The calculation of the first term at X+ is similar to the one in Eqs. (4.86)-(4.87): 

1 



T+[u,<t>,i — ►oo] f 



but the result is different. We obtain 



5/ 



(4.99) 
(4.100) 



(4.101) 



x+[«,</.,r^oo] 45 dv? 
with // in (4.80). Here the 1/r^ term survives but it has the same form as (4.81), and the 
same inference applies. 
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5 The early-time behaviours 



Let T be the support tube of the physical sources in (1.8). By the original assumption 
about the asymptotic stationarity of external fields, in the past and future of tube T (T~ 
and T"*" respectively) there exist asymptotic timelike Killing vectors such that all sources 
J in (1.8) are conserved along their integral curves. Let ^"(x) be a timehke vector field 
that interpolates between the KiUing vectors at T~ and T+, and up to 0[3?] is the Killing 
vector for the whole of T. Denote the Lie derivative in the direction of and 

J= (£^J'^^ C^P) . (5.1) 

• 

For simplicity, the support of J will be assumed compact. Then, on the central geodesic 
of any Bondi-Sachs frame, there will be two points o" and with u{o~) — u~, u{o'^) — 
and u~ such that the support of J is entirely inside the future light cone of 

o~ and the past light cone of o+. In the approximation (3.13) the fields of nonstationary 
sources propagate only between the cones u — u~ and u — rather than in the whole 
causal future of supp J. We have 

y^^i. ^ yz.^^ ^ Q ^ J^O, U<U-, U>U+ . (5.2) 

For what follows we need some properties that the functions involved in the calculation 
possess in presence of the KiUing vector field. 

Geometrical two-point functions in presence of the Killing 
vector. 

If is a Killing vector field, then the Lie derivative £^ commutes with the covariant 
derivative when acting on any object: 

£^y/3^a... ^ V^r^X"^- . (5.3) 
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Indeed, by definition, 

= ^'^V^X"- - EX'^-V^r (5.4) 
whence making use of the KiUing equation one obtains 

£^v^x"- - v^£^x"- = nx^- iy^v^c - ^"/.r) , (5.5) 

and this commutator vanishes by the known property of the KiUing vector field [25] 

VaV0i, = R^pJi. . (5.6) 
Consider now the world function a = (7{x, x) for a timelike or spacelike separation of 



the points x and x. The vector n — V(T/y2|cr| is a unit tangent at the point x to the 
geodesic connecting x and x, and the vector n — Va/^2\a\ is the oppositely directed unit 
tangent to the same geodesic at the point x. Since the quantity ^"'ria is conserved along 
the geodesic [25], we have 

Cna = -Cna , (5.7) 

or 

e"V„(7(x,x)+rV„a(x,x)=0 . (5.8) 

This is the conservation law for the world function. By continuity it holds also for a null 
separation of the points x and x. 

Eq. (5.8) can be written in the form 

(>C^ + £^-)(7 = (5.9) 

whence using (5.3) one obtains 

{C^ + C^)Va, . . . V„„V^, . . . Vp^a = (5.10) 

— a conservation law for the derivatives of the world function. 

A similar law for the propagator of the geodetic parallel transport can be obtained as 
follows. From the equations (3.3) for g""^ we have 



£5a^jV^^"" + a'^V^(/:eO = , (5.11) 
{Cia^)v,g''^ + a^V,{jC^g"^) = (5.12) 
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where use is made of the commutation law (5.3). Since 



by (5.10), combining Eqs. (5.11) and (5.12) yields 



(5.13) 



We have also 



+ Cc]g 



which follows from the initial condition in (3.3) and the condition 







(5.14) 



(5.15) 



(5.16) 



The latter condition is obtained by differentiating the equation for g""^ in (3.3) and setting 
X — X. With the initial condition (5.15), the solution of Eq. (5.14) is 



(5.17) 



This is the desired result. 

Finally, consider Eq. (5.8) with a timelike 4" and go over to the limit x in this 

equation. The insertion of the asymptotic behaviour of the world function (4.5) yields 

- eV^«(x) + ^'^V,r{x)^U^{x) - -e^,U^{x) , x^I+ . (5.18) 

Here the term with d/d(j) drops out by (4.9), and, with the usual normalization of the 
timelike Killing vector at infinity, we have 



1 . 



x+ 



(5.19) 



Thus we obtain 



e^M^) = 1 (5-20) 
which is the conservation law for the null hyperplanes. Using (5.3) we obtain also 



C^Va,...VaMx)^0 , n>l . 



(5.21) 
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In the Bondi-Sachs frame with the retarded time normahzed as in (1.14), the solution of 
Eq. (5.20) is 

U^{uJ,f)^u + L^{^,f) , (5.22) 
where the function L^{(j),f) possesses the properties 



d 

' a^^^^'^''^'' 







(5.23) 



by (4.13) and (4.14). 

It should be emphasized that the relations above for the two-point functions hold 
only in the case where the geodesic connecting the two points lies entirely in the Killing 
domain. 



Retarded kernels in presence of the Killing vector. 



That an arbitrary timelikc vector in Eq. (3.12) and similar equations has been denoted 
is no mere coincidence. As will be seen from the calculations below, it is advantageous 

to choose for this arbitrary vector the vector defined in the beginning of the present 

section. 

Consider an integral of some scalar source X over the past hyperboloid of a point x, 
and calculate 

j dxf^^6(a{x,x)-q)x (5.24) 

past of a; 

= / dxg'/'^^^C^j{aix, x) - q)x 

past of X 

= - / rf^^'/'5(<7(x,x)-g)v„(r||^x) . 

past of X 

Let now the observation point x belong to the past Killing domain u{x) < u" . Since 
the hyperboloid in (5.24) lies entirely in the causal past of x, the integration point x also 
belongs to the past Killing domain. Furthermore, the timelike geodesic connecting x and 
X lies inside the past light cone of x and, therefore, passes entirely through the Killing 
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domain. Then, using Eq. (5.8) and the corollary Va^"" = of the Killing equation, we 
obtain 

J dxf/H{a{x,x)-q)x= J dxf/^S(a{x,x)-q)C^^X , 

past of X past of x 

u{x) < u" . (5.25) 

For obtaining a similar result in the case of a tensor source X, consider X in Eq.(3.8). 
Write 

rV.X = C-^{g% . . . g^^^^)x-^---^-{x) + g>^l^ . . . g%C^X-^^-^-{x) (5.26) 
and use the law (5.17). It is then seen that Eq. (5.25) generalizes as follows: 

J dxg'/'S{a{x, x) - q)g%^ . . . = (5.27) 

past of X 

= I dx^V25(^(^^^)_^)^m ^ ^,(^)<^- . 

past of X 

Eq. (5.27) implies that, in the past Killing domain, the Lie derivative JC^ commutes 
with the operator Tig in (3.20): 

C^HqX{x) = rLqC^X{x) , u{x) < u- . (5.28) 

It then commutes with all the vertex operators: 

C^F{m,n)XiX2{x)^F{m,n)[C^Xi)x^^F{m,n)X^{C^X^ , (5.29) 

C^^F{m,n)X^X2{x) = ^F{m,n)(C^X^X2 + ^F{m,n)xAc^X2) , (5.30) 

C^-^F{m,n)X,X^{x) = -i— F(m, n) (^C^X^jx^ + -^F{m,n)xAc^X^ , 

I I2 I \\\ I2 ' '1' '2 

(5.31) 

u{x) < u~ 

and, since (5.28) is valid for g = as well, it commutes with the retarded operator l/D: 

C^^X{x) = ^C^X{x) , u{x) < u- . (5.32) 



56 



Causality of the vacuum radiation. 



For any of the physical sources J the commutation relations above yield 

C^^J{x) = ^Jix) = Q , (5.33) 
C{HqJ{x) = 7^, J(x) = , (5.34) 
C^F{m,n)JiJ2{x) = F{m,n) Ji J2{x) + F{m,n)Ji J2 (x) = , (5.35) 

u{x) < u~ 

and similarly with the other vertex operators. Here use is made of Eq. (5.2) and of 
the retardation of all kernels. The relations above and the commutativity of with 
the covariant derivative suffice to infer that the scalar I{x) in Eq. (2.37) possesses the 
property 

C^J{x) = , u{x) < u- . (5.36) 
Consider now the moment of the scalar /, and calculate 

^Dr(nA\I) = -^Jdxg^''S'{U4x)-u)i (5.37) 

On the integration hyperplane in (5.37) we have u — U^{x) > u{x) by (4.13). Therefore, 
with u < u" we have u{x) < u~ for the whole of the hyperplane. Furthermore, the null 
geodesic connecting the point x of the hyperplane with the point (f) at T"*" belongs to this 
hyperplane itself. Therefore, we may use Eq. (5.20) and the Killing equation to obtain 



u<u~ 



^ j dxf'H[u^{x) - u)C\/ J ■ (5.38) 



Hence, by (5.36), 



. (5.39) 



u<u~ 



Recalling that the flux of the vacuum energy is expressed through the moment in 
(5.39) by Eq. (4.20), we arrive at the following final result: 

dM{u) 



du 



= . (5.40) 

u<u~ 
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The classical radiation is also governed by the time derivatives of the moments but 
these are the moments of directly the physical sources J. Also for these moments we have 



= (5.41) 

u<u~ 



by (5.2). 

We proved the following assertions. Stationary sources radiate neither classically nor 
quantum-mechanically. Moreover, they don't produce even quantum noise. Radiation, 
including the uncertain oscillations of the energy flux, starts not earlier than the first light 
signal from a nonstationary source reaches the observer at infinity. 

The proof of these assertions given above is adjusted to the approximations made in 
the present paper but the assertions themselves are valid beyond these approximations. 
By using the commutation law (5.3) and the boundary condition of retardation, it is not 
difficult to show that, in the past Killing domain, the operator commutes with the 
exact retarded resolvent. Then it commutes with all retarded form factors of the form 

r dm? . . . dm^ p(m? , . . . m^) 
J (□i-mi^)...(n„-m„2) 

In other words, if the source and background fields in the equation (□ — ■m^)(fi — J are 
static in the past, then so is the retarded solution (f and so are all functions 

r(ni,...n„)Ji... J„ . (5.43) 

Furthermore, as remarked in Sec.2, expression (2.48) is exact, and / in this expression is 
generally a sum of terms (5.43). Owing to the presence of the overall time derivatives in 
(2.48), Eq. (5.40) is an exact fact. 

Convergence of the vertex operators. 



We may now come back to the question of convergence of the vertex operators in 
(3.23) and (3.32)-(3.33). The vertex operators act directly on the physical sources J and 
are expressed through the operator Tig in (3.20). By rewriting the integration measure in 
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(3.20) in the form 

dxfl' = (5.44) 



where dT, is the induced volume element on the hyperboloid a{x,x) = q , and using that 
on this hyperboloid (Vcr)^ = 2cr = 2g , we obtain 

HgJix) = - — ^== / rfSJ , g < . (5.45) 



47rv/=2^ 

past sheet of a{x,x)=q 

It is important that , with x fixed and q — oo, the entire past hyperboloid a{x, x) = q 
shifts to the past and finds itself in the Killing domain where the sources J are stationary. 
Since, in addition, the support of J is confined to a spacetime tube whose intersection 
with the hyperboloid a{x,x) = q remains compact as g ^ — oo , the integral in (5.45) 
tends to a finite limit. Therefore, with x fixed and q ^ —oo we have 



oc — ^ . (5.46) 

q—' — oo Y q 



It follows from (5.46) that, with the sources asymptotically static in the past, the 
vertex functions (3.23) converge for all m and n, and so do the functions (3.32) but the 
functions (3.33) diverge logarithmically for all m and n. This is a consequence of the 
singularity at the zero- mass limit of the operator !/(□ — m^)^ applied to static sources: 
the integral (3.26) with such sources diverges at = 0. 

The remedy is in the fact that the chain of retarded kernels connecting T^^^ with 
the sources J contains time derivatives. This is, in particular, the derivative d/du in Eq. 
(4.18). Up to higher-order terms in 3? it can be commuted with all kernels and considered 
as acting directly on one of the J 's in the guise of the Lie derivative C^. This is seen 
from Eqs. (5. 38), (5. 31), etc. which, up to higher-order terms in 3ft, hold everywhere by 
the definition of the vector field When dealing with the vertex function (3.33), this 
commutation should always be assumed done. 

If at least one of the two sources in the vertex function has the support properties 
of J, all vertex functions converge including the function (3.33). Indeed, the support 
of J is compact and, therefore, no matter where the point x is located 111 , at the limit 



^^If the point X is located at u{x) < u , all past hyperboloids of x are outside the support of J, Eq. 
(5.34). 
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q —oo the past hyperboloid of x will go out of this support (see Fig. 5). Let q(x) be 
the parameter q of the earliest hyperboloid of x that still crosses the support of J. Then 
Eq. (5.46) gets replaced with 



Tig J (x) 



= , (5.47) 

q<q(x) 



9 

and all integrals (3. 23), (3. 32), (3. 33) with one of the sources X — J acquire a cut off at the 
lower limit irrespectively of the nature of the other source. 

That the vertex function (3.33) with one of the sources X =J converges is seen also 
from the fact that any power of the masslcss retarded operator !/□ applied to J, (!/□") J, 
is well defined since, by (5.33), all functions (l/n"~^) J vanish identically at I~ . The 
integrand in (3.26) is then no more singular at = than in the convergent integral 
(3.25). 

The presence of J among the sources ensures the convergence of an arbitrary super- 
position of retarded kernels since any convergent integral of the form 

Y{x) = J dx^gl'^ . . . dXn9l^^G{x\xi, . . . x,) Ji X2 . . . X, (5.48) 

where G is a retarded kernel possesses the property 

Y{x) = , u{x) < u- . (5.49) 

Owing to (5.49), the function Y{x) vanishes identically not only at the past timehke and 
null infinities {i~ and I~) but also at spatial infinity (i^): 



Y(x] 



= . (5.50) 

x— >i~ or I~ or i'^ 



This is illustrated in Fig. 5 whence it is seen that, at any of these hmits, the past light 
cone of X will go out of the support of J . Owing to the property (5.50), any integral of 
Y{x) with a retarded kernel converges, and any convergent integral of the form 

Z{x) = j dxigl'"^ . . . dxngl/'^G{x\xi, . . . Xn)YiX2 ...Xn (5.51) 

possesses again the property (5.49): 

Z{x) = , u{x) < u- . (5.52) 
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It is worth emphasizing that the support of a function hke the Y[x) or Z[x) above is 
no more compact and, moreover, is not even a spacetime tube since it has a null boundary. 
This support (shown with broken lines in Fig.5) is the causal future of the support of J . 
Therefore, the presence of time derivatives in the kernels has generally no effect on their 
behaviours at the future infinities beginning with the future null infinity for u > u~. The 
integrals (4. 51), (4. 52), etc. are, however, cut from below owing to (5.41). 



61 



6 The late-time behaviours 



Eq. (5.40) proves that the total energy of vacuum radiation is determined indeed by 
the hmit of late time in (2.49). Prom (4.20) and (5.39) we have 

M(-oo) - M(oo) = - ^ <fS{<t>) D^M\I) (6.1) 

whence it follows that we are to study the behaviour of the moment Dx as m ^ oo, i.e. 
the behaviour of the retarded Green function (4.17) in the future of X+. For (6.1) to be 
finite and nonvanishing, this behaviour should be 

Dx{u,(t)\I) (xu^ . (6.2) 

The retarded Green function in the future of X+ . 

We shall show that the behaviour of 

Di(u,(l>\X) = Di(u,(l>) (6.3) 

as It — > oo is determined by competing behaviours of X{x) at the following four limits. 

i) The hmit of X{x) as u{x) — > oo along the timelike hues filling a spacetime tube 
(T). For such lines one can take the hues r = const., = const, of an arbitrarily chosen 
Bondi-Sachs frame. This limit will be denoted 

Xr+[r,(f),u ^ oo] . (6.4) 

ii) The hmit of X{x) as x moves to the future along the timehke geodesies that 
reach the asymptotically flat infinity. In the Bondi-Sachs coordinates these geodesies are 
asymptotically of the form 

" u — s , (f) — const. , s— >oo , 0<7<1 (6.5) 



where s is the proper time, and 7 is the boost parameter (1.22). This limit has been 
denoted in Introduction 

Xi+[y,(f),s ^ 00] . (6.6) 
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iii) The limit of X{x) as x moves to the future along the null geodesies 

u — const. , = const. , r — > oo . (6.7) 
This hmit has been denoted 

r — > oo] . (6.8) 

iiii) The limit of X{x) as x moves to the asymptotically flat infinity along the spacehke 
geodesies. In the Bondi-Sachs coordinates these geodesies are asymptotically of the form 

u = -(1 - f3)r , = const. , r^oo , -1</3<1 . (6.9) 

This limit will be denoted 

X,o[/3,0,r^oo] . (6.10) 

Consider the integration hyperplane in (4.19). One of its null generators is radial i.e. 
crosses the central geodesic of the Bondi-Sachs frame at some point o. As follows from 
the consideration in Sec.4 (and as illustrated by Fig.4), the integration hyperplane lies 
outside both sheets of the light cone of o. Therefore, the support of X in (4.19) is confined 
to the exterior of this cone. In Fig. 6 the hght cone of o is depicted with bold lines, and 
its exterior is divided into four subdomains. The support of X in (4.19) is their union 

suppX = lUllUlIlUlV , (6.11) 

and the integral (4.19) itself is a sum of the respective four contributions 

= d} + d}^ + + d}^ . (6.12) 

Subdomain 1 belongs to a tube r < tq with tq sufficiently large for subdomains 
II, III, IV to be already in the asymptotically fiat zone. Subdomain III is bounded by 
two future light cones u — Ui and u — U2 with large positive U2 — \u2\ and large negative 
Ui — —\u]\. With the Bondi-Sachs parametrization of the integrand in (4.19) we have 

ro oo 

Dl{u,(f)) = — j d^S{^) j dff'^ j du\{Vu,Vf)[^5{u^{u,^,f)-u) 



-oo 

oo oo 



dI^{u,<P) = ^ j d^S{<p) j dff'^ j du\{Vu,Vf)\ ^ 6{u^{u,<p,f) - u)x , (6.14) 
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dP^{u,(I)) = — J d^S{^) J dff^ J du\{Vu,Vf)\~^5(y^{u,^,f) -u)x , (6.15) 



OO wi 



^ j d^S{<p) Jdff^ J du\{Vu,Vf)\ ^6(u^{u,(p,f) -u)x .(6.16) 



The argument u of the function D\{u^ (j)) labels the future light cone of o. As -u ^ cxd, 
the point o moves along the central geodesic to the future. In addition, the parameters 
ro, wi, U2 should be made functions of u such that , as — > oo. 



ro{u) 



OO 





ui{u) 



-OO 



, 



U2[U) 
U2{u) 



+00 , 

. 



(6.17) 



u u u 

It is then seen from Fig. 6 that, as u — > oo, subdomain I shifts to T+, subdomain II shifts 
to i+, subdomain III shifts to and subdomain IV shifts to Let us show this by a 
direct calculation. 

Consider first the contribution of the tube, Eq. (6.13). At a however late u the 
hyperplane U^{u,(f),f) — u will cross the nonstationary region between the cones u — 
u" and u — , and, therefore, will not belong entirely to the future Killing domain. 
Nevertheless, at a sufficiently late u the intersection of the hyperplane with the tube will 
be at > It"*". The null generators of the hyperplane emanating from this intersection to 
the future will then also be at -u > -u^. Therefore, in the integral (6.13) with a sufficiently 
late u one may use Eq. (5.22) owing to which the integrand becomes restricted to 



u = u — Lij){4>,r) , M — > OO 



(6.18) 



Since, in addition, the range of the integration variables (p, r in (6.13) is compact, X turns 
out to be at T"*". We obtain 



1 1 

= — / d^S{(j)) [ dff^ (Vm, Vf) '^Xr+[f, 0, -u = -u ^ oo] , (6.19) 
3 47r J J 



and the geometrical factor in the measure is bounded by virtue of the original assumption 
that the metric has no horizons. 

In the contributions of the domains II, III, IV one can use the flat-spacetime expressions 
for the metric and the function U^{u,^,f), Eqs. (1.15) and (4.15). In the integral (6.14) 
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go over to the integration variable 7 according to the formula 

f = -^^u (6.20) 
1-7 

i.e. parametrize this integral with the family of geodesies (6.5). It follows that, when 
w — > 00 , 7 ranges in the interval < 7 < 1 , and, with 7 fixed, X turns out to be at 
i+. In the integral (6.15), u ranges in a bounded interval, and, with u fixed, X turns out 
to be at X"*". In the integral (6.16) go over to the integration variable /9 according to the 
formula 

f=-^. (6.21) 

i.e. parametrize this integral with the family of geodesies (6.9). Then /3 ranges in a 
bounded interval, and X turns out to be at In this way we obtain 



(6.22) 



1 — (1— COS (i))|u2 I/m 



A-K J J (1 — 7C0sa;)^ 



M — >• 00 



ro/u 



(1 — 7C0sa;)^ ' ' 1 — 7C0sa; 



1 [ (PS{(i)) f ,_/_2^\ r- 7 - 1 
/ V-ry I j„.l„2v\ ^u,(p,r=- ^00] 



(6.23) 



j du{f^x)_ 



Att J l — coscuJ ^ ^2:+ 1 — coso; 



(6.24) 



1 — (1— COS Lu)\ui\/u 



An J J {p — cosujV^ 'v^ p — cosu 



An J J — ' ' 3 — cosu 

COSOJ ^ ' 

Here u = uj{(f), 0) is the function of the angles defined in (4.16). 

By (6.17), the integration limits in (6.22) reach as u — > 00 the end points of the 
interval < 7 < 1. As explained in Sec.l, the end point 7 = 1 supports the contribution 
of the future of X+ and so does the end point 1*2(1*) — > 00 in (6.23). The behaviours of 
X at these end points are related by Eq. (1.29). Similarly, the end point 7 = supports 
the contribution of the boundary of the expanding tube T+ and so does the end point 
^o('^) ^ 00 in (6.19). Finally, the upper limit in (6.24) reaches as u — > 00 the end point 
(3—1 which supports the contribution of the past of T+. Another such contribution 
comes from the end point ui{u) — > —00 in (6.23). The end-point contributions are of 
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measure zero only if the respective integrals converge at these end points; otherwise, they 
may be essential. (The end point P — —1 which corresponds to the future of X~ makes 
no contribution.) 

For simplicity, we may confine ourselves to scalar sources X which are, moreover, 
stationary in the past: 

d 

—X(u,(j),r)=0, u<u- (6.25) 
au 

since the source / in (6.1) possesses these properties. The behaviour of X at is then 
restricted by the convergence condition (4.25). A restriction on the behaviour of X at 

can be obtained by considering the retarded Green function with a static source (the 
Coulomb potential). The convergence condition is in this case 



X 

Therefore, we may put 



X,o[/3,0,r^oo] = ^ , n>A (6.27) 

r , 1 f yl(0)/r" , U<U~ 

\B[u,(t)) /t'^ , u>u 

with some coefficients A and B. For X at and it suffices to consider the power 
behaviours 

Xr+[r,(j),u ^ oo] = C{r,(j))u'' , (6.29) 
X,+ [7,0,s^oo] - W{^,cl>)sP (6.30) 

with arbitrary k and p. However, to agree with (6.28) and (6.29), the coefficient of the 
latter behaviour should have appropriate singularities at 7 = 1 and 7 = 0: 

W(-f, 0) = (1 - 7)-^+^/^7^"^t«(7, (/)) (6.31) 

where w{'j, (p) is a regular function. By the correspondence (1.29) we have also 



= w{l,(l>)2P/^ uP+^ , (6.32) 



and similarly 

C{r, 0) 



^w{0,(Py-'' . (6.33) 
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Upon the insertion of the behaviours above in (6.19)-(6.24) and the use of (6.17) the 
following results are obtained. The contribution of i° is 

dP^{u,(P) =0{u-^) . (6.34) 

The contribution of proper, i.e. of a finite range of retarded time along T+, is also 

^0{u-^) . (6.35) 



The contribution of proper, i.e. of the open interval < 7 < 1 , is 



(6.36) 



Finally, the contribution of the tube proper, i.e. of a tube with compact spatial sections, 
is 

Dl{u,(f)) ^0{u^) . (6.37) 

As far as the end-point contributions are concerned, they become essential only in the 
case where the behaviour (6.36) is comparable with (6.35) or (6.37). If the power in (6.36) 
is different from the one in (6.35) and from the one in (6.37), then the total result is a 
sum of these competing powers: 



0(u-^) + 0(uP+^) + 0{u'') , p + 3^-1 , p + 3y^k, 



(6.38) 



and one is to pick out the dominant one. If, however, p-|-3 = — lorp-|-3 = /c, then the 

joint contribution of the coincident powers gets amplified by a factor of log u: 

n , , , / 0(ti-Mog ti) + 0(«') , p + 3 = -1 

''''"■^'.-oo=lo(„Mog„) + 0(«-),p + 3 = * ■ '-'■""^ 

A comparison of the behaviours obtained with the one in (6.2) shows that relevant 

contributions may come only from X at or T+. However, at , there arc three extra 

powers added to the exponent of the behaviour of X at late time, Eq. (6.36). If X 

doesn't grow at either of the limits (as is normally the case) or has one and the same 

power of growth at i'^ and T+, then only X at i'^ is capable of making a contribution to 

the real vacuum energy. Therefore, the most important case is generally the one where 

the contribution of X at is dominant: 



p> -A , p>k-Z 



(6.40) 
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In this case the result is 



M— »00 



— / d 5(0) / d-f- — — — Xi+[7,0,s = ^ oo] 



Air J J (1 — 7cosco')^ 1 — 7COSC1; 

(6.41) 

with coscj in (4.16). With the specification in (6.31) this integral converges. 
Formula for the energy of vacuum radiation. 



It will be shown below that for X = I the conditions (6.40) are fulfilled with p = k, 
and the exponent p in (6.30) has the needed value p = —1: 

J,+ [7,0,s^oo] = iiy(7,0) . (6.42) 

Then Eq. (6.41) yields 
Di{u,(P\I) 



1 2 



^=—J d'S{^) I d7-^=L=^(l-7cosa;(0,0))' iy(7,0) . (6.43) 

According to (6.1), the coefficient of in this expression is the angle distribution of the 
vacuum radiation. 

The integration over the directions of radiation, i.e. over the angles at 0, can 
be done explicitly: 

I d'S{<P) (1-7 cos ^(0, 0)) = JY^^ (6-44) 
whence for the total radiation energy (6.1) one obtains the result 

M(-oo) - M(oo) = / ^77^ / d'Si<l>) .^^I'.t, (6.45) 

\ I ) 

in terms of the coefficient in (6.42). Eq. (6.31) for the case (6.42) is of the form 

l^(7,0) = (l-7)'/'^(7,0) (6.46) 



^^Not to be confused with the angles </> at i+ on which the function W{'y, (f) in (6.42) depends. A 
summary table of the integrals over the 2-sphere used in the paper is given in Appendix B. 
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which ensures the convergence of the integral (6.45). The coefficient W{'y,(j)) in (6.42) 
will be calculated and shown to be negative definite. 

It will be recalled that condition (6.46) emerges from the correspondence between the 
hmits i+ and and is a consequence of the convergence condition (4.26) at T+. However, 
this correspondence concerns only the future of X"*". Therefore, for the validity of (6.46) 
and hence for the validity of the result (6.45) it is only important that the convergence 
condition at hold at late time. There should be a value u' of retarded time, no mattar 
how late, such that 

(/ ) =o(^) . (6.47) 

The vertices and trees in the future of X+ . 



The analysis above of the behaviour of Di{u,4>\X) as w — > oo will now be applied 
to the case X — J where J is any of the physical sources (1.8). Since the support of 
J is confined to a tube, of the four contributions (6.34)-(6.37) there remains only the 
contribution of T"*", Eq. (6.37). Since, at T"*" , J is stationary, we obtain immediately 



Di{u,cf>\J) 



(6.48) 



This result can be made more precise. Let u' be the value of retarded time such that 
for u > u' the intersection of the hyperplane U^{u,(t)^f) = u with the support tube of 
J is entirely in the future Killing domain u > . Then for u > u' one may repeat the 
calculation in (5.37)-(5.38) to obtain by (5.2) 



|-Di(«,0|J) 



. 



u>u' 



(6.49) 



Eq. (6.49) determines the behaviours of the vertex functions in the future of X"*" since 
it provides a cut off for the time integrals. Thus from (4.29) we obtain 



{F{m,n)JiJ2 



X+[u,(/),r— »oo] 
oo 



(6.50) 



2(-l) 



m+n+l 



U 



m+n 
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Specifically, the behaviour of the vertex function V^caiar in (2.43) is 



(v; 



scalar 



) / 

i/u— >oo J 



dr 



(6.51) 



X+K</),r-»oo]^'"-*°° r _^ 

where the abbreviation (4.36) is used. In the same way, for t4ross we obtain from (4.64) 

lu 



:/9 
cross 



X+ [■«,</), r—>c»] 



dr 



§-2Di\r,<P)][^,IDf{r,4>)]. (6.52) 



More important are, however, the behaviours of the vector and tensor vertices since 
they are 0(l/r^) at . For V^ect we have from (4.65) 



(VaV;3K:; 



:/3 
vect 



\ 1 U 

]/ ti— >oo 6 



rl rl 1 

/ dry„4-I)f(r,</.)][-^/(r,0)]+O(-) (6.53) 



but the result for Krav is completely different: 



grav 



(6.54) 



owing to (4.81) and (6.49). The latter result applies also to the only tree that is 0(l/r^) 
atX+, Eq. (4.101). 

Inspecting Eqs. (2.37)-(2.39) one infers that the behaviours obtained here and in Sec. 4 
can be summarized as the following result for the scalar /: 

(i.M.^=o,).>.' ^ 6^^ / '^^^ - + 0{^) (6.55) 

where the only 1/r^ term at late time comes from the vector source. Since, in (6.55), 
T < u and 



tr 



< 



(6.56) 



^dr ^^dr 

(see Sec. 8), the integrand is negative definite. The only possibility for vanishing of the 
1/r^ term in (6.55) is, therefore. 



tr 



= 



'dr ^ ^dr 

By imposing this limitation on the vector source we force condition (6.47) to hold. 



(6.57) 
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The trees at and . 



The inference from the analysis above is that the real vacuum energy is determined 
by the behaviour of the scalar I at i'^ and possibly at T+. We begin the study of these 
behaviours with the trees. 

The only nonlocal object that figures in the trees is (!/□) J(a;) with !/□ in (3.13) and 
J{x) in (1.8). The past light cone of x with x moving to or T"*" is similar to the null 
hyperplane shifted to the future in a sense that there exists a value u' of retarded time 
such that for u{x) > u' the intersection of the past light cone of x with the support tube 
of J is entirely in the future Killing domain. Then also the null geodesies connecting the 
points of this intersection with the vertex x of the cone are entirely in the Killing domain. 
Therefore, at u{x) > u' one may use Eqs. (5.32) and (5.2) to obtain 

£^lj(a;) = , u{x) > u' . (6.58) 

For the behaviour at this gives straight away 

0{u°) . (6.59) 



□ 



T+ [r,</),w— >c»] 

As far as the behaviour at i'^ is concerned, it is the same as the behaviour at X+ for 
u{x) > u' since, in the static case, there is no dependence on the direction at infinity. In 
this way we obtain 

-J =0(-) . (6.60) 

n i+[7,</,,s-»oo] 

It follows that at r+ all the trees (4.82), (4.84) and (4.96) are 0{u^). At i+, the 
trees (4.82) vanish identically, and the trees (4.84) and (4.96) are 0(l/s^). In terms of 
Eqs. (6.29) and (6.30) the maximum exponents of the trees are, therefore. A; = and 
p — —2 whereas for making a nonvanishing contribution to the real vacuum energy these 
exponents should beA;>2,p>— 1. Thus the contribution of the trees vanishes. 

Having excluded the trees we may go over to the chief thing: the behaviours of the 
vertex functions at 
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7 The late-time behaviours (continued) 
Spacelike hyperplanes. 



The behaviour of the world function as one of its points tends to i"*" and the other one 



stays in a compact domain can be obtained as follows. Since —2(t{x, x) is the geodetic 
distance between x and x, in the principal approximation as x — > i"*" there should be 



yj-2a{x,x) ^s{x)(l + 0) , O 



(7.1) 



where s{x) — > oo is the proper time of the point x moving to i^. In a curved spacetime 
there is one more growing term, proportional to logs(a;). Therefore, generally we have an 
expansion of the form 

^-2a{x,x) =c(.)-4^ + 0(-) , (7.2) 

c{s)^s + 0{\ogs) (7.3) 

where c(s) is restricted by an additional condition that it doesn't depend on x, and a 
notation is introduced for the 0{s^) term of the expansion. In this term, T^^{x) is some 
function of the point x depending also on the coordinates 7, of the point x, and the 
normalization factor \/l — 7^ is introduced for further convenience. The T^^{x) is defined 
up to an addition of an arbitrary finite function of x i.e. a function of 7, 0. 

The insertion of the asymptotic behaviour (7.2) in the equation (3.2) for a with respect 
to the point x yields 

(vr(x))' = -(l-7') , T{x)^T,^{x) , (7.4) 

and, since Va{x^ x) is past directed, so is S/T{x). The function T^^[x) with fixed 7 and (p 
defines the family of spacelike hypersurfaces 

T^^{x) ~ T — const. (7.5) 

with time r growing towards the future. The vector field orthogonal to these hypersurfaces 

Ar"(x) ^ y^ Ifd^ , N\x) = -l (7.6) 
VI - 7 
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is a gradient and at the same time has a unit norm. A combination of these two properties 
signifies that the integral curves of this vector field 



ds 



(7.7) 



are geodesies. The specification of these geodesies can be read off from Eq.(7.2): 

- Ar"(x) = V°^-2(t(x,x) . (7.8) 

[7,(}!),s— »oo] 

At every point x the vector field N"'{x) is tangent to a timelike geodesic that, when traced 
to the future, appears at the asymptotically fiat infinity with the energy E = {1 —7^)^^''^ 
at the point (p of the celestial sphere. In a compact domain the geodesies having one 
and the same 7 and (p differ by "translations" and make a 3-parameter congruence. It 
follows from (7.2) that the geodesic congruence thus defined is hypersurface-orthogonal, 
and the orthogonal hypersurfaces are just (7.5). The hypersurfaces (7.5) will be called 
hyperplanes, and there are different families of spacelike hyperplanes for different values 
of 7 and 0.0 

By using Eq. (6.5), the expansion (7.2) can be rewritten in terms of r{s) 00 : 



—2a{x, x] 



[7,(/),s— »oo] 7 

whence, going over to the limit 7 — > 1, we obtain 



r + 0(logr^ 



(7.9) 



cr X, X 



7-^1 



i4> 



+0 



7=1 



o 







(7.10) 



By the correspondence (1.29), this sequence of limits should coincide with the future of 
On the other hand, at 1^ we have the expansion (4.5) in which the function U^{x) 
is independent of the time along X+. It follows that 



UJx) + const. 



7=1 



(7.11) 



where the const, is independent of x. Thus we infer that the function T^^{x) admits 
the limit 7^1 and, at this limit, the spacelike hyperplanes (7.5) turn into the null 



^^The Bondi and ADM masses refer to the Lorentz frame at infinity in which the asymptotically flat 
spacetime rests as a whole, i.e. the center-of-mass frame. The parameter 7 refers to the same frame. 
Therefore, the dependence on 7 cannot be boosted away and is real. 
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hyperplanes (4.7). The geodesies orthogonal to the spacehke hyperplanes also become 
null and turn into the generators of the null hyperplanes. 
The insertion of expansion (7.2) in (3.4) yields 

V7 -^7(?^('^) 



5/(a;,a;)V^s(x) 



x-*i+[j,(p] 



(7.12) 



which is the law of parallel transport of the tangent vector along a timelike geodesic 
orthogonal to the hyperplanes. By using the asymptotic form (6.5) of the geodesic, one 



may expand Vs over the vector basis at i"*" : 

1 



V^s(a:) = ^==^(v^i(x)-7V^r(x)) , x^i+[^,cf)\ 



Then the law of parallel transport takes the form 



(7.13) 



(7.14) 



and at 7 = 1 goes over into (4.11). 

In presence of a timelike Killing vector the spacelike hyperplanes possess the properties 
similar to (5.20) and (5.21): 



>C^-Vai • • • Va„T^</,(x) = , n > 1 



(7.15) 

(7.16) 



Finally, in the case where x is at the future asymptotically flat infinity or i+, one 
can use the flat-spacetime formula for T^^{x): 



x—*T+,i 



u + r(l — 7cosa;(0, 0)j 



^ t- -fni{(j))x' 



(7.17) 



with the same notation as in (4.15). 



The radiation moments and conserved charges. 



The function T^^{x) may formally be regarded as a two-point function with one point 



at 



T^^{x)=T{x,x) 
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(7.18) 



With the aid of this notation and with the sources in (1.8), define the following radiation 
moments: 



b{x) 



= ^ y dxf/^S(T{x, x) - t)p{x) 
^ dxf/'S{T{x,x)-T)g%ix,x)fi{x 
= ^1 dxg'/'5{T{x,x)-r)g°'^{x,x)g''^{x,x)J^~^{x) 



(7.19) 
(7.20) 
■ (7.21) 



These objects are respectively a matrix, a matrix vector, and a tensor defined, however, 
only for x at i'^. The specification i+ will, therefore, be omitted and we shall write simply 
D, i)", D""^. The similarly defined moments for the Ricci scalar R and the matrix Q will 
be denoted and D'^. In addition to the parameters 7 and (j) of the point at i+, 
the radiation moments depend on the parameter r of the hyperplane. As follows from 
(7.11), the previously introduced moments Di are the special cases of the moments D 
corresponding to 7 = 1: 



7=1 



7=1 



7=1 



, etc. 



(7.22) 



Owing to the law of parallel transport (7.14) and the conservation laws (1.9), the 
vector and tensor moments satisfy the relations 



(Vai-7Var)^I)" = , 
(Vat-7V«r)^L'"^ = 



(7.23) 
(7.24) 



which generalize (4.42) and (4.48), and are proved similarly. The conserved quantities in 
these relations 



e = {V^t - -/W^r)b'' = dxt'H{T{x,x) - t)Wj,TJ~^{x) 



(7.25) 



and 



P"" = ^(V^t - lV,r)D^'' = ^J dxg'/'5{Tix, x) - t)/,(x, x)V^TJ^^(x) 



(7.26) 



are the full charges (4.44) and (4.50) since the latter are independent of the choice of the 
hypersurface E = crossing the support tube of J. 
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By using (3.6) and the conservation laws, the vector and tensor moments can be 
brought to the forms 



with 



dr ' dr^ 



^= -1-J dxg^/'^S(T{x,x) -T)a"{x,x)Wp,TJf'{x) 



V"f^{x) = — I dxf/^5(T{x,x)-T)a''{x,x)a^{x,x)VnTVpTJ^^{x) 



{1.11) 

(7.28) 
(7.29) 



Although we called moments the D 's, this name is more appropriate for the V 's. Their 
relation to the textbook multipole moments is considered in Sec. 8. The V 's always 
appear differentiated with respect to r. Undifferentiated, they are infinite at the limit %^ 
as distinct from the D 's. The consideration below is carried out in terms of the D 's in 
(7.19)-(7.21). 



The vertex operators at i+ . 



An essential point concerning the vertex functions is that the hypcrboloid o{x. x) = q 
with q ranging from to — oo sweeps the whole causal past oi x. As x ^ i^, the 
hyperboloid of x with any fixed q will go out of the support of the nonstationary sources 
but there will always be a range of q for which the hyperboloid crosses this support. For 
q to stay in this range, it should be shifted to — oo simultaneously with x i+. This 
suggests the following replacement of the integration variable in (3.23), (3.32) and (3.33): 



-2q = c{s) - -== (7.30) 
VI - 7 

where r is the new integration variable, c(s) is the function in (7.3), and 7, s are the 
parameters of a; G i~^['y,(j),s — >• 00]. With q replaced as in (7.30) and r fixed, we obtain 
from (7.2) 

5(a{x,x)-q) -1 s(T^^{x)-r) . (7.31) 



Hence 



nqX{x) 



^/^^7^ 

-D^{tA\X) (7.32) 
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where 



D^(r,0|X) = ^/ dxg'/^ 6 (r^^ix) - t)x 



(7.33) 



is the full (7- dependent) radiation moment of the test source X. The moment in (4.19) 
is its limiting case 



Di{T,cj>\X) = D^{T,cf>\X) . (7.34) 

7=1 

Another essential point concerns the spacelike hyperplane (7.5) at late time r. This 
case is analogous to the case of the null hyperplane. As the hyperplane (7.5) shifts to the 
future, its intersection with the support tube of J turns out to be in the future Killing 
domain as well as the timelike geodesies emanating from this intersection towards i'^. Eq. 
(7.15) can then be used to obtain for the spacehke hyperplane the result analogous to 
(6.49). Namely, there exists an instant r' of time r such that 



^D^(t,0|J) 



. 



t>t' 



(7.35) 



One is now ready for obtaining the vertex functions at i'^. Making the replacement 
(7.30) in the integral (3.23), and using (7.32) and (7.35) one arrives at the following result: 



m+n— 3 



F{m,n)JiJ2 

1 \ m+n—l 



(7.36) 



m+n+3 



dr 



/ d \ m+1 , I ^ T / d \ n+l ,i^s 

(-) L>^(r,0|JO][M D,{t,4>\J2) 



This is the sought for growth in time. Comparing this result with Eq. (6.50) one can 
check the fulfillment of the correspondence (1.29) between the limits and X"*". A similar 
calculation for the operators (3.32) and (3.33) yields 

^-F{m,n)JiJ2 = 

I i+[7,(/),s— ♦00] 



(7.37) 



I 0/ l\m+n— 1/ / \ m+n+1 P r / d \ m+1 ir / d \ '«— 1 



d \ n-l 



1 
2 



F{m,n)JiJ2 



i+[7,(/>,s— »oo] 



(7.38) 



„ / 1 \ m+n— 1 / / \ m+n—l f r / u \ m— 1 , , , n r / U \ n— 1 , , ■ 
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Specializing to the vertex functions in (2.40)-(2.43) one obtains 

— oo 
— oo 

— OO 

— oo 

In each case, only the highest exponents m and n of the vertex operators F{m,n) work. 

As seen from the expressions above, the function V^scaiar has precisely the needed power 
of growth at (cf. Eq. (6.42)) but the remaining functions grow apparently too fast. On 
the other hand, in Eq. (2.39) these latter functions appear differentiated. Unlike in the 
case of X+, here the derivatives alone help as we show below. However, the conservation 
laws simplify the result greately. 

The non-scalar vertices at z+ . 



scalar 



i+[7,0,s-^oo] 



yal3 
' cross 



0/3 



vect 



i+[7,0,s^oo] 



j+ [7,<)!>,s— >oo] 



In view of (7.13), the conservation laws (7.23) and (7.24) for the moments can be 
written in the form 

V«s-^D" = , V„s-^D"^ = . (7.43) 
dr dr 

Therefore, when projected on Vs, the senior asymptotic terms of the vertex functions 
(7.40)-(7.42) vanish. For further use we shall record this fact in the form 



where V""' is any of the vertex functions bearing indices 



0(-\/"-) (7.44) 



""^^When referring to tensors at infinity we always mean their components in the nuU-tetrad basis or, 
equivalently, in the Minkowski frame. The components of radiation moments in this frame are finite. 
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For tensors V"''^ and V^f^^^ possessing the property (7.44) we are to calculate the 
behaviours at i"*" of the quantities 



(7.45) 



appearing in (2.39). For this purpose it is convenient to go over to new coordinates of the 
point at i+: 

X = (s, 7, 4>) = (s, f ) , y = in\(t>) (7.46) 
i+ 

where s,7* are the new coordinates, and n*(0) is defined in (1.27). The coordinates s,7* 
can be related also to the Minkowski coordinates at . In the Minkowski frame (1.25) 
the geodesic (6.5) is of the form 



t 



s — >• oo 



(7.47) 



which gives the transformation law. 

The transformation law (7.47) makes it easy to obtain the metric at in the new 
frame. We have 



= (Vs, Vs) = -1 , g^' ^ (Vs, V7^) = , 



(7.48) 



where Qik is the 3-dimensional metric depending only on 7': 

1 /c , lilk \ r k 



7 ^ - 7 

-ife /'I ^!2\ts:ik ^i^,k\ 



(7.49) 



r-(l-y)(^^ -7^) , V^^Vdet^: 



1 



/ifc 



(1 - 7^)2 

Covariant derivatives with respect to the 3- metric g^k will be denoted Vj , and the following 
relation will be quoted: 

^i^k—7=^ = ^ Qik ■ (7.50) 
VI - 7 V 1 - 7 

It will also be noted that, in terms of 7^ the integral over the parameters at i+ that 
figures in (6.45) is of the form 

1 



^77^ j (fS{(t))X 



«+ [7,(/),s->oo] 



i+[7',s— ►oo] 



:7.5i) 
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By a direct calculation in the metric (7.48) one obtains 

v„v,i-<-'' = v,v.y*+(4+4)fcy»+(2|+2)v.v«+(^+5|+l)v'». (7.52) 

which is valid identically for any symmetric V^^. That every index "0" on V " is accom- 
panied here by an extra power of 1/s (or d/ds) is a general rule. A similar calculation 
for a symmetrized fourth-rank tensor yields 



ds 



+ terms in -V^^mn J_yOOmn J_yOOOn 



1 



0000 



(7.53) 



S S"" 

In Eqs. (7.52) and (7.53) there figure the coordinate components of V^- i.e. the 

projections 

yo- = VaS 1/"- , 1/^- = Vaf 1/"- . (7.54) 

The virtue of the introduced coordinate frame is in the fact that the set V^s , Vq7* is not 
a finite vector basis at as distinct from the coordinate basis of the Minkowski frame. 
From the transformation law (7.47) it follows that 



0(1) , V„7' 



(7.55) 



Therefore, the projection of V""- in (7.54) that has n indices "i" has also the extra 

as compared to V"' - in the Minkowski frame. With regard for this fact, all terms of an 

expression like (7.52) or (7.53) are, generally, of one and the same order at i"*". We have 



VaVpV^^ -0{-V) , Vo.VpV,V,V^^^\ -Oi-V) . (7.56) 

i+ 5 j+ S 

This is precisely what is needed for the vertex functions (7.40)-(7.42) since inserted in 
(7.56) they all become 0{l/s). Specifically, the spatial components of all vertex functions 
in the frame (7.54) are 



ik 



0(1 



yik 
1 vect 



ikmn 
grav 



0(1 



(7.57) 



The inference above obtains without using the conservation equation (7.44). Using 
this equation one can assert more, namely that all terms in (7.52) and (7.53) with at least 
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one index "0" on V'" can be discarded since they decrease by one power of 1/s faster than 
the "i" terms. Using also (7.57), one obtains for both V^ross K^ect the expression 



(7.58) 



and for V^rfv expression 
V«V^V^V,l^("^'^^) = (v,VfeV„V„ + 26^,feV^V„ + 45^,fe5mn)^^*'™"^ + 0(4) (7-59) 

i+ S 

These expressions are next subject to the integration over the parameters at i"*" with 
the measure in (6.45): 



1 

/ d,,^ J d^Si<P)^^-^^i. ..)=! d^.^g Oil - 7) ^^(. . .) . 



(7.60) 



Therefore, one needs them only up to terms vanishing in this integral. Relation (7.50) 
makes it easy to integrate by parts in (7.60) or, equivalently, to bring expressions (7.58) 
and (7.59) to the following form: 



afi 



(7.61) 
(7.62) 



with some quantities and in the total derivative terms. It can be checked that these 
quantities vanish indeed at 7 = 1 which is the boundary of integration in (7.60). 

One may now come back to the covariant form. With the metric in (7.48) one has 

s%kV'^- = {gap + V„sV^s)l^"^- , (7.63) 

and the projection VaS^" " vanishes at i'^ by virtue of the conservation equation (7.44). 
Thus we obtain the following final results: 



al3 



-^dapV^^ + total derivatives , 
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(7.64) 



VaV^V^V.l^^"^'^") = -TQapg^uV^"^''''^ + total derivatives . (7.65) 



Substituting for and \/"^^^ the expressions (7.40)-(7.42) we obtain the behaviours 
at i"*" of all non-scalar vertices (the total derivatives are omitted): 

= -7^(1-7^^ / dr(-^ 

i+[7,</),s-»oo] 



1 °° /72 j2 
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—oo 

i~^\^.d).s-^oo\ 1ZU6 J Z 



^dr^ )^dT^ 



With the behaviour of V^scaiar in (7.39) included, we have the complete hst. All the be- 
haviours are of the form (6.42) i.e. proportional to 1/s, and for all but one the proportion- 
ality coefficient contains the factor (1 — 7^)^/^ as required in Eq. (6.46). The exception 
is the vector vertex (7.67) in which instead of (1 — 7^)^/^ we have (1 — 7^)^/^. This is a 
consequence of the lack of convergence at X+. That the results are correct can be checked 
by comparing the asymptotic expressions above with the behaviours in the future of T"^, 
Eqs. (6.51)-(6.54). The correspondence (1.29) between the hmits i'^ and X"*" holds in all 
cases. 

An alternative way of handling the non-scalar vertices is considered in Appendix B. 
The vertices at . 



Since we assume that all timelike geodesies are infinitely extendable to the future, 
the behaviour (7.1) holds at T+ as well. The parameter \J —1q is then to be shifted by 
s like in Eq. (7.30), and, as a result, there will appear integrals over the hypersurfaces 
orthogonal to the geodesies coming to T+. These integrals analogous to the moments will 
be stationary at late time like in Eq. (7.35) since their only ingredients will be the metric 
and the stationary source. It follows that, in terms of the proper time along the world 
lines filling the tube, the behaviours of the vertex operators at are the same as at i^: 



F(m,n)JiJ2 



Since we assume that there is no future horizon, the proper time of an observer in the 
tube is analytically related to the external time u (normalized in (1.14)). Hence 



F{m,n)JiJ2 



ock"^+'*-^ , u^oo . (7.70) 
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Eq. (7.70) signifies that the vertex terms in the scalar I have one and the same power 
of growth at and whereas, for making a contribution to the real vacuum energy, 
the growth at T"*" should be by three powers of u faster (Sec. 6). Thus we infer that the 
only contribution to the total vacuum energy comes from the vertex functions at i^. 



^^For 14calar this follows directly from (7.70). For the non-scalar vertices the overall derivatives in the 
expressions like VcV^gV^g^!^ , etc. have at T+ the same effect as at This is clear from the treatment 
of these derivatives in Appendix B. 
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8 Creation of particles and radiation of waves 
The energy of the vacuum particle production. 



The result of consideration in the previous section is that formula (6.45) can indeed 
be used, and 1^(7, 0) in this formula is the sum of expressions (7. 39), (7.66), (7. 67) and 
(7.68) with the factors 1/s detached: 



W^(7,0) 



1'^ 

1 



2\5/2 



dr tr 



d 



d^ 



^dr^' 



1/ rf2 



+ 



3(1 

i 



^dr 
d' 



30 



By (1.12) 



b - -D^ 
6 



^dT^ 



D 



(8.1) 



(8.2) 



where D is the radiation moment (7.19) of the potential matrix P. One may notice that 
the coefficients of the first two terms in (8.1) are precisely such that the cross contribution 
PxR cancels. This cancellation (although nontrivial) is a consequence of adding — to 
the potential in Eq. (1.2) and is the only manifestation of the conformal properties of the 
effective action [31]. Otherwise, these properties play no role in the present calculation 
as the reader can see. 
Using also that 

= -g^^.Df^" (8.3) 
(cf. Eq. (4.37)) and inserting (8.1) in (6.45) we obtain the following final result for the 



^"There is an unceasing controversy in connection with the old paper [32] about the significance of 
the trace anomaly in dimensions higher than two. The calculation in the present paper may serve as a 
commentary to the following conversation that took place between two persons: 
N: The anomaly is a window through which we can see ... 
V: ... the back yard. 
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quantity (1.37) in terms of the radiation moments (7.19)-(7.21): 

1 oo 

^77^ [ dr [ (fS{(j)) tr 



M(-oo) - M(oo) 



1 



{Any 



D 



1 



i.4) 



There is a pole at 7 
hmitation (6.57): 



1 in the term with the vector moment. However, under the 



tr 







7=1 



.5) 



'dr ^ ^dr 

the pole cancels. As shown below, this limitation is a condition that the vector connection 
field contains no outgoing wave. 

Expression (8.4) is the total energy of the particles produced from the vacuum by 
external fields. It is useful to realize that the integral over 7 in (8.4) is none other than 
the integral over the energies of the outgoing particles. The vacuum radiation is thus 
obtained along with its spectrum. 



Positivity. 



Owing to the conservation laws for the moments, all contractions with the metric in 
(8.4) are positive definite. In particular, by (1.23) and (7.24), 



1 



(8.6) 



+ 2(1-7^ 



|^D"^V„r + \{^,D-'rn^m*, - 2(1 - 7^)|^D"^V.rV,r 



which proves the positivity of the gravitational-field contribution in (8.4). 

The positivity of the matrix contributions follows from assuming the self-adjointness 
of the operator of small disturbances of the quantum field . There should exist a 



^-'^For simplicity, let this field be boson and non-gauge. In the general case the calculation will include 
diagonalizing and squaring operators, and combining loops [24]. Eq. (8.4) gives the contribution of a 
generic loop. 
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symmetric and nondegenerate matrix ujab such that the operator 

ooAcH^ (8.7) 

with in (1.1) is the Hessian of an action. The self-adjointness then requires that i) 
the operator (8.7) be symmetric, and ii) the matrix cuab be positive definite. The latter 
is a condition of the absence of ghosts. The symmetry of the operator (8.7) imphes that 
i) i^AB behaves hke a metric with respect to the covariant differentiation: 

V^a;AB = , (8.8) 

and ii) ljab converts the potential matrix into a symmetric form: 

P^cjcb-PE^ca^O . (8.9) 

From (8.8) it follows that 

[V^, V.]u;ab = -n%,u;cB - n%^,ucA = , (8.10) 

i.e. loab converts the commutator curvature into an antisymmetric form. Using (8.8) once 
again one obtains the antisymmetry relation for the source of the commutator curvature 
in (1.7): 

J^"o;cB + Jb'^^ca = , 4" = > . (8.11) 

At this point it is important that the integrals with matrices contain the propagators 
of parallel transport for the matrix indices. The exphcit forms of Eqs. (7.19) and (7.20) 



are 



jAa 

'b 



(8.12) 
(8.13) 



Di = ^Jdx-g'/'6(T{x,x)-T)g\g/P, 
^b" = ^Jdx-g'/'6{T{x,x)-r)g\g/g%J 
By (8.8), the law of parallel transport for loab is 

coABg\g''B=^AB ■ (8.14) 

Using this law it is not difficult to prove that the symmetries of the sources in (8.12) and 
(8.13) imply the symmetries of the moments: 

D^cucB^D^cucA, (8.15) 

D^-u;cB^-D^''u;cA. (8.16) 
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The latter symmetries remain unchanged under a differentiation with respect to r and a 
contraction of (8.16) with any vector. 

Let F"*" and F" be any matrices possessing the properties 



r^^ ujcB ^ ±r^E ^cA . (8.17) 

The positive-definite metric uoab can be expanded over an orthonormal basis: 

ooAB = hA{M)hB{N)6{M,N) , hA{M)h'^{N) = 6{M, N) , /i^(M) = a;-^^^/iB(M) 

(8.18) 

with S{M,N) the Kronecker symbol. Denoting 

F±(M, N) = h^{M)T^^hB{N) (8.19) 

one has 

tr(f±)' = F±iF±^ = ±F±ia;^cr±^a;-^^^ = ± ^ (f±(M, TV))' . (8.20) 

N,M 

Thus, identifying r='= with 

where ko, is an arbitrary vector, we obtain 

' • (8-22) 

The first of these inequalities proves the positivity of the potential contribution in (8.4). 
Finally, by (1.23) and (7.23), 

!l^>(ib-)(lb^) = {io''n:){lD^,ni)+{l-,%lb'v4 . (8.23) 



Hence, by the second inequality in (8.22), 

which proves the positivity of the commutator-curvature contribution in (8.4). 

The matrix ujab itself does not figure in any of the final expressions (8.22)-(8.24). 
Only its existence is important. 
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Radiation of waves. 



Generally, the flux of the non-coherent radiation caused by pair creation in the vacuum 
is only one term of the mass- loss formula (1.35). To clear up the meaning of condition 
(8.5) and to complete the discussion of the radiation moments, we shall consider also the 
other terms which in the field-theoretic case stand for a radiation of waves ( not necessarily 
classical [23]). 

By analogy with the electromagnetic field we may assume that the contribution of the 
vector connection field to T^^^^ in (1.32) is of the form il 

^source = " ^tr (^„^7^^"7^''^ - ^^^'^7^"^7^„^) + a contribution of the potential P (8.25) 

where the change of the overall sign as compared to the case where TZ'^'^ is literally replaced 
by the Maxwell tensor is owing to the antisymmetry (8.10). The density of the flux of 



^source through X+ is then 



1+ 



47r ^4 
1 

47r 



(8.26) 



m. 



1+ 



2 -P--- ■■-yv2 

where the last form is obtained by inserting expression (1.21) for Qajs and using the 
antisymmetry of 7^"^. If we define the complex news function H for the waves of the 
vector connection field 



d 

— Cvect(M,0) 

ou 



1+ [?i,(A,r— >ool 



(8.27) 



then, according to (8.26) and (1.35), the quantity 
is the energy flux of the outgoing radiation of this field 



d 2 

— Cvect(M,0) 
OU 



(8.28) 



Outside the support of the source of TZ^^" . 
^"^This term was introduced in [26,27] for the gravitational waves but the electromagnetic waves can 

be discussed along the same lines. Two real components of the news function arc the initial data for the 

two independent degrees of freedom of a radiation field counted per point of X. 
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Assuming that there is no incoming radiation, we may use the Jacobi identities to 
express TZ"-^ through its source as in Eq. (A. 17) of Appendix A. Using also (4.17) we 
obtain to lowest order in 3? 



V"u ^Df{T, 0) - Vf^u ^D^{r, 0) 
ar ar 



Hence 



and 



d d - 

— Cvect(T, (t>) = mo, — D^{T, (j)) , 



d 2 

— Cvect(T, 0) 



d " ^ 2 

dr 



(8.29) 



(8.30) 



(8.31) 



where, in the last form, are the spatial components of D-[' in the Minkowski frame, 

and use is made of Eq. (1.28). Using the basis decomposition (1-21) for the metric and the 
conservation law (4.42) for the moment, expression (8.31) can be written in the covariant 
form 

|^C.ect(r,0)f=,.,(^^i")(|:^i'^) . (8.32) 
The total energy of the waves of the vector field emitted for the whole history is, therefore. 



— OO 

It is now seen that condition (8.5) can be written in the form 

tr 



(8.33) 



d 2 

— Cvect(T,0) = , 



(8.34) 



and its meaning is that the source of the vector field radiates no waves. 

Only the 7 = 1 moments Di appear in classical radiation theory as is clear from Eq. 
(4.17). However, even in this theory the concept of boosted (7 - dependent) moment D is 
useful since it helps to understand the nature of the multipole expansion. In the case of a 
nonrelativistic source, 7 is the only parameter that contains the velocity of light since 7 is 
the velocity of a particle at per unit velocity of light. Therefore, in the nonrelativistic 
approximation, D at 7 = 1 may be calculated as an expansion at 7 = 0: 



7=0 07 



+ ... 



7=0 



^.35) 
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This expansion gives rise to the multipole moments. Since at 7 = the dependence on the 
direction at should disappear, the expansion in 7 is also an expansion in the direction 
vector rf. For the spatial components of in the Minkowski frame one may write quite 
generally 



7=0 



7=0 



dr 



(8.36) 



■1 rf2 



d 



^A-r-T^ d^''^ + T- d^''^) + a term oc n' , rf^^^^fe = (8.37) 
^ dr^ dr ' 



and so on. Here d\ d^'''\ d^'''^ are the dipole, quadrupole and magnetic moments, and the 
term oc ra* drops out of (8.31). The appearance of the time derivatives in expressions 
(8.36) and (8.37) is owing to the fact that the multipole moments are coefficients in the 
expansion of the V in Eq. (7.28) rather than the D. For the radiation of a nonrelativistic 
source one obtains the standard textbook result 



1 

47r 



d^Si^) 



d 

— ( 

dr 



^VCCt V 5 



2/ 



3Wr2 



d'] + 



d^ 



180 



2 1 / rf' 



3Vc/r2 



d^^^^) +... . 

(8.381 



The (differentiated) news function for the outgoing gravitational waves is 
^Cgrav(M,0) = -\Va,vV^vmpm,(rR'^^'^ 



1+ [u,(f>,r^oo\ 



(8.39) 



Assuming again that there are no incoming waves and solving the Bianchi identities to 
express i?"^'^'^ through its source [15,23] we obtain to lowest order in 3? 

^2 



{rR 



X+[t,4>,t~*oo\ 



.40) 



where both types of brackets [] and <> denote antisymmetrization of the respective 
indices. Hence 



^;^Cgi.av(T, • 



^m^mp^Dfir,^) . (8.41) 



With the initial condition (5.41) one may integrate this equation to obtain the news 
function 

(8.42) 



r^^Cgi.av(''", > 



1 d T^aPi 



*See, e.g., [23] and the original references [26,27]. 
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Thus 



1 

" 2 



2 _ 1 

~ 4 



2 



(8.43) 



((^im - ninm){Skn " JT'fe'T'n) " ^((^jfe " nink){Smn " JT'mJT'n) ( J-^i*") (^--Or") 



where, in the last form, are the spatial components of D^^ in the Minkowski frame. 
Bringing this expression to the covariant form yields the result 

|-C,rav(r,0)[ = \{9^,gp. - l9a,9,.){i^Df){j-Dr) ■ (8.44) 

Therefore, according to Eq. (1.35), the total energy of the gravitational waves emitted 
for the whole history is 



ijdrj d^«s(0)l(,.,,,. - 1,.,,,.) (Id'^^) (Id 

— oo 

For a nonrelativistic source. 



7=1 



(8.45) 



7=0 



+0(7) , 



and 



D 



ik 



7=0 



(T'^ + a term oc 5 



ik 



ik 



(8.46) 



(8.47) 



where d^^ is the quadrupole moment, and the term oc 5*^ drops out of (8.43). Here again, 
0?**^ is the coefficient in the expansion of the "D in Eq. (7.29) rather than the D. Therefore, 
in (8.47) there appears the second time derivative. For the radiation of the gravitational 
waves by a nonrelativistic source Eqs. (8.43) and (8.47) yield the textbook result 



47r 



d^S{(t)) 



Ac 



(t,0) 



45Vdr3 



+ ... 



(8.48) 



Compare Eqs. (8.33) and (8.45) with (8.4). The similarity is striking. As a result of 
the derivation above, the quantum problem of particle creation becomes almost the same 
thing as the classical problem of radiation of waves. The whole difference is that, in the 
quantum problem, there figure the boosted moments D and, instead of setting 7 = 1, one 
is to integrate over 7. In the nonrelativistic case even this difference disappears since all 
moments will be expanded at 7 = , and the integral over 7 will be removed trivially. 



91 



Note also that in the case of the vector field both the classical and quantum radiation 
effects are determined by the second-order time derivative of the moment (in terms of the 
V , Eq. (7.27)). In the case of the gravitational field there appears one more distinction 
between the classical and quantum contributions to the radiation energy. Namely, the 
classical radiation is determined by the third-order and quantum by the fourth-order time 
derivative of the respective moment T> . This is a consequence of the dimension of the 
coupling constant in the quantum loop. 
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9 Specializations and examples 
Spherical symmetry. 



The fact that the energy of vacuum radiation is expressed through the moments with 
7 different from unity has an important consequence. Inspecting Eqs. (8.6), (8.23) and 
(8.43) one can see that at 7 = 1 there survive only the projections of the moments on the 
2-sphere S (the transverse projections) whereas at 7 ^ 1 there are also projections on Vr 
i.e. on the direction of motion of the partice at i"*" (the longitudinal projections). For a 
spherically symmetric source, the vector moment can have no projection on the sphere, 
and the projection of the tensor moment on the sphere can only be proportional to the 
metric on the sphere in which case the orthogonality (m, m) = comes into effect. It fol- 
lows that spherically symmetric sources cannot emit waves but can produce particles from 
the vacuum. This makes spherical symmetry an interesting case for studying the effect of 
vacuum radiation. In addition, the limitation (8.5) is in this case fulfilled automatically 
and so we can study the particle production by vector fields. 

For spherically symmetric sources we have from (8.23) and (8.6) 



'dr ^ ^dr ' ^dr 

- ls,.S.,){^D'^) i^D^") = li^D" + 3(1 - j')^D'^V.rV,r 

(9.2) 

Since these quantities are spherically symmetric scalars, they do not depend on the angles 
at i'^. Therefore, the integral over (f) in (8.4) removes trivially, and we obtain 



M(-oo)-M(oo) = ^/^777c^rtr[(A_i,)^_l(i_I,"V.r 

-00 

+ + 3(1 - ^)^^D'-^V„rV,rf] (9.3) 

where there is no more pole in 7 including in the term with the vector moment. 
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The longitudinal projections of the moments are conveniently obtained as follows. One 
differentiates the law of parallel transport (7.14) with respect to 7. This yields 



^i^r{x) g''p_{x, x) 



^ -^M^^70(^) + i^A^) - 7V,r{x))^g^^{x,x) 



a;6i+[7,0] 

(9.4) 



The derivative 

i:{9%{x.x) )=0[^ (9.5) 

can be calculated as a result of parallel transport of a vector along the closed contour 
consisting of two geodesies which emanate from x and come to with boosts 7 and 
7 + 57. With the contribution of (9.5) neglected , one obtains from (7.20) and (7.21) 

V«rD" = / dxtlH{T^^{x) - t)(v^|-7;^)>(x) + 0\^^\ , 



The metric of a spherically symmetric spacetime is generally of the form 



(9.6) 
(9.7) 



dS'' , dV' = gAB{y)dy^dy^ , dS^ = gab{<l>)drd<f>' 



La J lb 



(9.8) 



where F is some 2-dimensional Lorentzian space, S is the unit 2-sphere, and r = r{y) is 
a function on F. The line r = is a boundary in F. We shall denote g^ the determinant 
of gAB{y) and the derivative in F. The world function on a spherically symmetric 
spacetime is a regular function of the following arguments: 



a{x,x) — a{y,y, cos cu) , uj — uj{(t),(f)) 



(9.9) 



where \uj \(f), (j)) is the world function on the unit 2-sphere, and an expression for cos a; is 



given in (4.16). Hence 



T^4>{x) = T-y{y, cos a;) , a; = a;(0, 0) 



and Eq. (7.4) takes the form 



1 



d 



Denote 



g^^'WA^^B^ + ^(1 - cos^u;)(- 

H ^ a cos a; 



T± = 7;(y,±l) 



=-(1-7^) 



(9.10) 



(9.11) 



(9.12) 
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By (9.11), the integral curves of each of the vector fields 

N^\y) = -^V%- , N_\y) = -^V%- (9.13) 

are timelike geodesies in F. With fixed 7, the geodesies generated by A'"+ make a 1- 
paramctcr congruence foliating T . The geodesies generated by make another such 
congruence. The geodesies from the two congruences may be joined pairwise, a pair 
consisting of the two geodesies that hit one and the same point of the boundary r = 0. 
A world line in F consisting of two geodesies joined at r = is the mapping on F of a 
4-dimensional radial geodesic (the broken line in Fig. 7). This line is generated by iV+ in 
the future from r = and by in the past from r = 0. 

The spacelike lines orthogonal to the geodesies generated by also make a 1- 
parameter congruence foliating F. The spacehke lines orthogonal to the geodesies gener- 
ated by N_ make one more such congruence. These congruences are respectively 

T/^T , T-^T (9.14) 

with T the parameter of the congruence. By (9.11), 

d 



dcosou 



(TJy, cos u) ) = . (9.15) 

^ r{y)=0'' 



It follows that the spacelike lines (9.14) join pairwise at r = as shown in Fig.7. The 
lines Z!y+ = r and — r with one and the same r hit one and the same point of the 
boundary. The hyperplane Tj^{x) — r maps on the interior domain between these two 
lines. However, the boundary of this domain, i.e. the line 

(T/=r)u(T^-=r) (9.16) 

itself corresponds to no 4-dimensional original. Only in flat spacetime is this line a map- 
ping of a spacehke geodesic. 

In two limiting cases, 7 = 1 and 7 = 0, there come about coincidences. At 7 = 1 , the 
line = r coincides with a geodesic generated by iV+ , and the hne = r coincides 
with a geodesic generated by N_ . Both geodesies are then null, and their union is the 
mapping on F of a 4-dimensional radial hght ray. At 7 = , the line = r coincides 
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with the hne = r . In this case the mapping of the hyperplane on T degenerates into 
a hne. 

When all external fields are spherically symmetric, there appear only moments of the 
sources X that either do not depend on the angles or depend on them through uj{(t),4>) 
where are the angles at i'^. The moment of any such source can be calculated as follows: 

D^^j dxt'H(T^^{x) -T)x = -j (fy^f^f^ j d{cosuj)6(T^{y, cosu) - t)x . 

(9.17) 

Here the equation of the integration hyperplane is to be solved with respect to cos a;. Let 
cos a; = / be the solution. Then the moment is obtained as the integral in F 







-1 

X 


dcOSLU 







(9.18) 



cos UJ=f 

in which the integration domain is bounded by the lines T^" — r and T^'^ — r. This 
calculation will be carried out below for fiat spacetime, and here it will only be noted that 
the Jacobian appearing in (9.18) is 



dT^ 



dcosu 



smo; 



(9.19) 



where L{y,u!,^) is the angular momentum of the geodesic that, when traced backwards 

in time, comes from i+[7,0] to the point (y, 0) of a compact domain. This follows from 
the geodesic equation (7.7) whose angle component takes the form 

du! lid 



ds- r^v^W^-^-^^"'"'"^ ■ 
Electrically charged shell expanding in the self field. 



(9.20) 



It is not the purpose of the present paper to discuss applications of the result (8.4) 
but to give a simple example we shall consider the particle production by an electrically 
charged spherical shell expanding in the self field. Below, e, m, and S are respectively the 
shell's charge, rest mass, and energy in excess of the rest energy. 
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In terms of time t orthogonal to r , (Vt, Vr) = , a spherically symmetric electro- 
magnetic field is completely determined by a single arbitrary function e(r, t) which is the 
charge contained on the hypersurface t — const, inside the 2-sphere of area 47rr^: 

e{r,t)^ j dxg^/^S{t-t)9{r -f)Vpif , (9.21) 

V^F'^^ = 47rj^ . (9.22) 

Here F^^^ is the Maxwell tensor, and is the electromagnetic current. The choice of 
e(r, t) is limited only by the condition of regularity of the field at r = , e(0, t) = , and 
the condition that the support of belongs to a spacetime tube. Let r = p{t) be the 
equation of the boundary of the tube. Then 



e(r, t) 



— e — const. (9.23) 

r>p{t) 



where e is the full charge. In fiat spacetime, with the normalization of time (Vt)^ = — 1, 
the general solution of the conservation equation is 

- -(V^O^^ - (V-^)^^ . (9.24) 

and the general solution of the Maxwell equations is 

F""" ^ E^Vt - E'^Vn , E^^^^^V^r (9.25) 

where only an electric field is present. 

In the case where the support of the charge is a thin shell, we have 

e{r,t) =ee{r - p{t)) (9.26) 

and, by (9.24), 

M^,t)f = ^^S{r-p{t)) , (9.27) 

where r = p{t) is a law of motion of the shell. In this case there is only a static Coulomb 
field outside the shell; inside the field vanishes. However, this simplicity is only apparent. 
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Because the shell should expand in the self field, this field is nonstationary, and the law 
of expansion is not very simple il: 

mc 1 e 9 ^ , , 

mc^ + S . (9.29) 



The shell can be contracted as much as one wants by raising its energy but for given 
energy its minimum radius is 

rmin = ^ . (9.30) 

Set free, it expands monotonically from r = rmin to r = oo with increasing speed. 
We shall consider two limiting cases, the nonrelativistic shell 



dt V P 
and the ultrarelativistic shell 



^ = J-(2S--) , S«mc^ (9.31) 



P = ^(l-^) + yc2(t-t3tart)2+ (^^!^^)^ , t3tart = const. , £»mc^. (9.32) 

It is important that in the latter case the shell will remain timelike. The timelike ultra- 
relativistic shell is always slow near r = rmin and only at r ^ oo it approaches the speed 
of light. 

Particle production by a spherically symmetric electromagnetic 
field. 



The expanding spherical shell emits no electromagnetic waves but, as we show below, 
it excites the vacuum of charged fields, and the vacuum emits the quanta of these fields. 

The coupling of an external electromagnetic field to the charged vacuum fields will 
be introduced through the following form of the covariant derivative in (1.2): 

= d^i + qA^n , trCl^ <0 (9.33) 

^^The appearance of 1/2 in the Coulomb energy is owing to the fact that the force exerted on the 
surface charge is determined by one half of the sum of the electric fields on both sides of the surface [33] . 
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where q is the charge of the vacuum particles, is a numerical matrix, and the negativity 
of the trace of is a corollary of the self-adjointness of the equation of the quantum field 
(see Sec.8). 

For the commutator curvature and its source, Eq. (9.33) yields the expressions 

7^'^^ = g^F'^^ , P^A-Kqhf (9.34) 

in terms of the Maxwell tensor and the electromagnetic current. It is convenient to factor 
the coupling parameters out of the longitudinal projection of the vector moment: 

VarZ)" = qnD\\ . (9.35) 

Then for D\\ we obtain from (9.6) 

D\\^- j dxg^lH{T^^{x) - r) {v -^—%^)f . (9.36) 
With the metric of the Lorentzian section in (9.8) 

dV^ = -de + dr'^ (9.37) 
one may use expression (7.17) for the hyperplane: 

T^(f,{x) = i — fj cos u!{4>,4>) . (9.38) 

Hence 

r^± = t>^f, (9.39) 

and the hnes (9.14) are spacehke geodesies (Fig. 7). With a spherically symmetric current 
we obtain 

(Vp-^f^^)F = - cosu{Vf,f)r , (9.40) 
and Eq. (9.36) takes the form 



oo oo 



D\\ = J d^S{4)) J dff^ J dt 6{t - T - f 'J cos uj) cos uj{V fif)f . (9.41) 

-oo 

The integration over the angles in this expression brings the longitudinal moment to its 
final form 

oo oo 

D\^ = ^ J dtj dr e{t-T + r^)e{T -t + r^){t- T){y ^r)f (9.42) 



' -oo 
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which is a speciaUzation of (9.18). Here the longitudinal projection of a spherically sym- 
metric current is to be inserted from (9.24): 

In terms of D\\ the vacuum energy released by a spherically symmetric electromagnetic 
field is, by (9.3) and (9.35), 

2 1 oo 

M(-oo) - M(oo) = ^(-tr n') I I dT[—D\)'' . (9.44) 

-oo 

Radiation of the nonrelativistic shell. 



The longitudinal moment of the charged spherical shell is obtained by inserting ex- 
pression (9.28) in (9.42): 

Here i+ and i_ are solutions of the following equations: 

t+-^p{t+)^T , i_ + 7p(i_) = T. (9.46) 

These are the time instants at which the world line of the shell intersects the lines (9.14). 
The moment (9.45) is an integral along the world line of the shell between the points t_ 
and t+ as shown in Fig. 7. Integrating by parts in (9.45), and using Eq. (9.46) we obtain 

For the nonrelativistic shell, the moment can be calculated by expanding it at 7 = 0. 
Since, by (9.46), 

i+ = T + 0(7) , i_=T + 0(7) , (9.48) 

this amounts to expanding p{t) at t = r. A straightforward calculation using only Eq. 
(9.46) yields 
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whence 



67 (P 



/(r) + 0(f) . 



(9.50) 



Inserting this expression in (9.44) one obtains the result 



M(-oo) - M(oo) 



2 2 

q e 



— ^3 

i-trCl') J dr(^p2(r))' (9.51) 



180 X 127r 

— oo 

which holds for the nonrelativistic shell irrespectively of the specific form of the law p(t). 

It is now important to know the entire history of the shell. We shall assume that 
before some time instant t — igtart the shell was kept at the point of maximum contraction 
f = ^min and next was let go. Beginning with t = tstart it was expanding unboundedly in 
the self field. This world line is shown in Fig. 7. 

Since 

dn(t) 

= , 



dt 



*<*start 



we have 



dt 



2\lll 



^start 



dp 

7 



iP 



2\l/l 



(9.52) 



(9.53) 



where the primes denote the derivatives of p{t) , rmin is given in (9.30), and the replacement 
of the integration variable t pit) is to be made with the law p{t) in (9.31). It is 
convenient to write 



(p^)'" = 6pV' + 2pp'A(p") 



dp 



and use the law of motion (9.31): 



P 



32 1 



\2E 



One obtains 



and 



dp 

7 



2m p- 

2\W 



2 ' P 



(P 



V m mp 



2 ryy^ 5 /2 



dx- 



X 



9/2 



/ 



dx 



X 



9/2 



2 2 2 
2 

3 5 7 



(9.54) 



(9.55) 



(9.56) 



(9.57) 



Finally, restoring h and c, we infer that for the whole time of expansion the nonrela- 
tivistic shell looses the energy 



M(-oo) - M{oo) = 



(-trQ2) g2 2^ 5/2 



81 X 25 X Tn hc^mc^ 



-(- 



'28 , 8 « mc^ . 



(9.58) 
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Radiation of the ultrarelativistic shell. 



Coming back to the exact expression (9.47) for the longitudinal moment, and using 
(9.46) we obtain 

rli . 1 rli 1 

(9.59) 



1-7 



dp 
'dt 



dt- 



and 



d 



e 
2^ 



1 



1 + 7 



dp 
'dt 

1 



where the notation assumes that the substitution t — t+ or t — t- is to be made. 

Since the law p{t) is different for t < 4tart and t > 4tart , the range of integration over 
T should be divided into three intervals, 



(9.60) 



D 



d 2 

dr + Jdr + l dT){—D\\^ 

— OO T\ T2 



(9.61) 



according to the location of the points t+ and t- : 



-OO < T < Ti 
Ti < T < T2 

T2 < T < OO 



t— ^start ) t+ ^start 
t— ^start ) ^+ ^ ^start 
t— ^ ^start ) t+ ^ ^start 



The motion of the points t+ and i_ along the world line of the shell as r increases can be 
traces in Fig. 7. We have 



drilo.f = , 



-OO 

T2 



1 2 



Kr 
and, from (9.46), 



d ^2 
Du 



474 J 



1 



'dt' 



(9.62) 
(9.63) 

(9.64) 



"T"! — 4tart ~ T'^min , 7'2 — ^start + T'^n 



(9.65) 
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In (9.63) and (9.64) the law p(t) pertains already to the expansion stage. 

Consider a shell that moves all the time with the speed of light. For such a shell, 
assuming that it passes through r = rmin at t = tstart , the solutions of Eqs. (9.46) are 

/ , N T — tstart + '^min dp{t) /n ««^ 

P{t±) = 7— , —77- = 1 . 9.66 

1 =F 7 clt 

and expression (9.60) vanishes identically. Hence we infer that a null shell creates no real 
energy from the vacuum. §1 Now consider a shell that before t = tgtart is at rest and then 
suddenly starts expanding with the speed of light. This shell creates an infinite amount 
of energy owing to the jump of its velocity at t = tgtart • Indeed, in this case the term 
(9.64) vanishes 

ldr{-D,y^Q (9.67) 

T2 

but the term (9.63) doesn't: 

n 

The latter expression behaves like 0(1/7^) at 7 = which in view of the presence of 7^ 
in the measure in (9.44) is a regular behaviour. However, at 7 = 1 expression (9.68) has 
a pole, and the total energy (9.44) diverges: 



(m(-oo) - M(oo)) oc -£ log(l - 7) 



(9.69) 

7=1 



The pole comes from the lower limit in (9.68) at which the regime changes from static to 
null. These considerations suggest that, in the case of the timelike ultrarelativistic shell, 
the effect may come only from a neighbourhood of r = rmin and it should be finite owing 
to the continuity of the shell's velocity. It is clear in advance that the divergent term 
(9.69) will be regularized with the regularization parameter mc^ / £ , and the regularized 
term will be of order £ \og{£ /mc^) . Since this term grows as mc^ /£ — , it will make a 
dominant contribution to the energy of the particle production. 
Indeed, for the timelike shell, Eq. (9.67) gets replaced with 

^dr(^DX^^ , ^-0 (9.70) 



T2 



^^This is not the case for a gravitationally charged sheU contracting in the self field. 
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while the term (9.63) can be expressed directly through the law of motion (9.32) by 
making the change r — >■ i+(T) of the integration variable and inserting the Jacobian from 
(9.59): 



1 1 



1 2 



/ 



(9.71) 



2 P{t+ij2)) o /I 

[ dp^ \ 2 ^(l-7p') 



2 

min 



J p'Lp2(l-7p') prjnin 

^min 

Here the approximation p' — 1 can be used in all terms except in the denominator (1 — 7p') 
where the deceleration of the shell near r = Tmin cannot be neglected. Specifically, the 
upper limit in (9.71) can be calculated in the approximation p' — 1 using (9.66) and 
(9.65): 

p{U{r2)) = T^r^in , ^ - . (9.72) 
i — 7 o 

However, for the insertion in the denominator (1 — 7p') , the law of motion (9.32) should 
be approximated better: 

p = , ^ = ■ (9.73) 

Y (P - Tmin)^ + (rmin-^)^ 



In this way, introducing the integration variable 

X = ^^^^ , (9.74) 

'"min 

we obtain for mc^/£ — > 

I dri^oX = -^A— ^ + 21ogi^ + 27 (9.75) 

J ^dr 'l^* 27ni-7 1 + 7 + ^ ' ^ ' 



7 /mc^\2 ' r ' , 1 



27/(1-7) 



/ dx— 
J 1 

^ 



I-7V £ ^ J (1 + x)2[(l - 72)x2 + (mcV£) 



21 



This expression has already a regular behaviour at both 7 = and 7 = 1 (cf. Eq. (9.68)). 

Upon doing the integral in (9.75) and next the integral over 7 one is to retain the 
senior terms at the limit mc^/S . These are seen to be the senior terms of (9.75) at 
the limit 

(1 - 7) - , ^ - , (1 - l){:^y = finite . (9.76) 
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At this limit expression (9.75) becomes 

oo oo r. 

Hence, interchanging the integrations over x and 7 , we obtain 

1 00 00 2 

jiW I - ^(.og^) /d. + 0(1) . !^ ^ (9.78) 

-00 ^ ' 

where the contribution comes from the end point 7 = 1. 

Thus we infer that for the whole time of expansion the ultrarelativistic shell looses the 
energy 

M(-oo) -M(oo) = ^4^^|^flog^)f , £»mc' . (9.79) 

This is already interesting since, at a sufficiently large £ , the shell radiates more energy 
than it has. Accounting for the backreaction of radiation becomes necessary. 

Both quantities (9.58) and (9.79) are independent of the charge e of the shell and 
remain nonvanishing at the limit e = . This paradox is a consequence of the property 
of classical theory that the charge cannot be switched off gradually. The backreaction 
of the vacuum may change this result since the expectation-value equations contain a 
dimensionless parameter e/q . Below we show that the range e/g^l is within the limits 
of validity of the approximations made. 

Validity of the approximations. 



The calculational scheme presented uses two types of expansion. The effective action 

is expanded in the number of loops, and within each loop order it is expanded over the 

basis of nonlocal invariants like in Eq. (2.2) i3 . The limits of validity of the latter 

expansion are presently to be considered. These limits are best seen from a comparison 

of the technique of nonlocal form factors with the Schwinger-DeWitt technique [30] . 

^^In the axiomatic approach of Refs. [19-22] the expansion in loops is avoided, and expansion (2.2) is 
considered as an ansatz for the full action. 
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As shown in [34], the nonlocal form factors emerge as a result of a partial summation 
of the local Schwinger-DeWitt series for the effective action. At a given order in the 
dimensional coupling constants this series can symbolically be written down as follows: 

S,^ = const. / dxg'/' ^ 4^ E(V)2^-2"(3?)" (9.80) 

p A* n=l 

where /i is the mass of the vacuum particles. The nonlocal expansion of the effective action 
is obtained by summing all terms in (9.80) with a given power of 3? and any number of 
derivatives. Within each order in //^ one neglects 3?"+^ as compared to . Therefore, 

the condition of vahdity of the resultant nonlocal expansion is generally of the form [34,19] 

3?^ « VV3? . (9.81) 

The partially summed series admits the limit — > [34,14-17] and at this hmit takes 
the form (2.2). There is no problem of principle in generalizing the nonlocal expansion to 
the case of massive vacuum particles but the neglect of the mass is correct if 

A(^3? « VV3? . (9.82) 

The approximation in which Eqs. (9.81) and (9.82) hold may be quahfied as the high- 
frequency approximation meaning the frequency of the external ( or mean ) field. 

Not all limitations implied in (9.81) and (9.82) should necessarily be fulfilled since not 
all terms neglected in the partial summation of the series (9.80) may be important for 
the problem in question. Thus, the effect of particle creation is caused by the external 
field's variability in time, and its spatial inhomogeneity is unimportant. Let I be the 
characteristic spatial scale of the source of the external field and u be its characteristic 
frequency. Then, for the problem of particle creation, the important limitation stemming 
from (9.82) is 

fjLC^ «hu . (9.83) 

Under this condition the results of the technique of nonlocal form factors are valid for 
massive vacuum particles as well. 

When specializing to the problem of particle creation by an external electric field one 
is to consider only the mixed components of the commutator curvature TZ/^i, . For these 
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components condition (9.81) yields 



noi«j , (9.84) 



and, by (9.34), 



Uoi = const. j-E (9.85) 



where E is the characteristic strength of the external electric field: 

E=^ . (9.86) 

Thus, for this problem, the condition of validity of the nonlocal expansion is 

^ « hi^ (9.87) 

t 

where the quantity on the left-hand side is the Coulomb energy of the vacuum particle 
in the external electric field. Condition (9.87) is to be combined with the condition of 
validity of the one-loop approximation: 

I- « 1 . (9.88) 
nc 

For the nonrelativistic shell, I and u are unrelated: 



e 



2 



g3/2 



^ = T ' ^ = ^7^ ' £«mc' . (9.89) 

In this case the energy of the shell is clutched by the inequalities (9.87) and (9.89) in the 

interval 



^-« /— 

he q V ?7T.C^ 

Owing to condition (9.88) this still leaves room. It is only important that the ratio e/q 
be not too large: 

2 

. (9.91) 

he q 

The weaker the coupling of the quantum field, the bigger charge of the external field can 
be considered. Thus, for the electron-positron vacuum, e = lOg is admissible but e = Vilq 
is not. 

For the ultrarelativistic sources we have 1/1 = 0, and the inequality (9.87) takes again 
the form (9.91). Thus condition (9.91) is universal but the relativistic case is better suited 
for the approximations made since no upper bound on the energy £ emerges. 
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Concluding remarks 



This section will not be long since it may be called " reminiscences of the future" with 
the exception of the following point. 

As remarked in Sec. 2, Eqs. (2.38)-(2.45) pertain to the case where 7^^/? and P are 
metric independent but the final result (8.4) is valid for arbitrary local TZap and P . An 
inspection of the table of asymptotic behaviours in [18] shows that, in the contributions 
of the variational derivatives STZap/^g^y and 5P/5g^i, to T^^^ , the highest exponents 
m and n of the vertex operators F{m,n) are by one smaller than in (2.40)-(2.43). The 
contributions of these variational derivatives are, therefore, pure quantum noise. This 
remarkable fact gives the result (8.4) the status of a generating expression applicable to 
fields of any spin [24] . 

The remaining limitations imposed on the external fields require bigger efforts for their 
removal. 

The limitation (8.5) signals that the theory contains another effect: the vacuum screen- 
ing or amplification of the electromagnetic waves emitted by a source. This is equivalent 
to an observable renormalization of all multipole moments. The emergence of the lim- 
itation on the vector field is connected with the fact that, in the case of this field, the 
nonlocal vacuum polarization and the local charge renormalization are of one and the 
same order by dimension. The separation of the respective terms in the effective action is 
quite subtle (this issue is discussed for QED in [35]). Hence it is clear why no hmitation 
like (8.5) emerges in the case of the gravitational field. It is also clear what should be 
done. Since there appear nonlocal sources that behave at X+ hke 0(l/r^) rather than 
0(l/r^) , two amendments are needed. First, terms 0(n°) of the form factors can no 
more be discarded. They should be extracted from the exact form factors [17] and added 
to (2.36). Second, expression (3.14) for the kernel of log(— □) is no more valid and should 
be replaced by the exact expression [21,22]. It can be expected that, after introducing 
these amendments, the difficulty with the convergence at I'^ will be removed along with 
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the limitation (8.5), and the result will be the following expression for the density of the 
radiation flux from the vector source: 

(cf. Eqs. (8.4) and (8.33)). Here Z is the missing renormalization constant. This effect 
is analogous to the effect of the vacuum gravitational waves [23] . The difference is only 
in the ways in which these effects appear in the formalism and in the dimensions of the 
coupling constants. The dimension of the coupling constant causes that the gravitational 
effect does not boil down to a mere renormalization [23]. 

Removing the limitation that the metric contains no horizon is a separate problem 
since in this case the total radiation energy is infinite but the reason why it becomes 
infinite can be pointed out straight away. What becomes wrong is the assumption of 
asymptotic stationarity of the sources in the future which leads to Eq. (7.35) and hence 
to Eq. (7.36). Since r is an external time, as r — oo the sources moving in the tube hit 
the event horizon. Therefore, the asymptotic stationarity gets replaced by the condition 
that at r ^ oo the sources remain finite along with their internal derivatives. The growth 
of the vertex function at i+ should then become by one power of s faster. 

Finally, the assumption that the support of the physical sources is confined to a space- 
time tube may seem technical and minor but is in fact physical and inadmissibly restric- 
tive. As pointed out in Sec.4, the significance of this assumption is in the fact that it 
excludes the radiation of charge. In the external-field problem this covers interesting cases 
but in the self-consistent problem this is merely wrong since what has been calculated 
above is exactly the radiation of the gravitational charge. Obtaining the backreaction 
equations is impossible without letting the sources out of the tube. 
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Appendix A. Identities for the vertex operators. Use 

of the Jacobi identities. 



The functions A^^{D^, □„) and B^'^{n^, □„) in the asymptotic expansions of the form 
factors (2.33) are given in Ref.[18] as hnear combinations with rational coefficients of the 
vertex operators (2.44) 



p., d ^kr^f d ^^fenlog(j„,□„,/j„□„) 



Expression (A.l) can also be written in the form 



, rriy^n . (A.l) 



F(k k)-nk^nk„j d Nfc^/ d xfc„ log(D^/DO 

The results in [18] simplify drastically if one takes into account the constraints that 
exist between the functions (A.l) with different km and kn . The constraint equation 
reads 

F{kra, K + l)+ F{km + 1, ^n) = -{km + K + l)F{km, k^) (A.3) 

and is obtained by acting with the operator 

ntnti^ti^t (A.4) 

on the easily verifiable identity 

(□ g g )Mn../nn)^_log(n./n„) ( 

The first step in simplifying the functions Aj^{Um., □„) and B^^{Um., □„) is a removal 
of the dimensionless coefficients Dnj/Dn that some of these functions have in front of 
F{km,kn) [18]. This removal can be done so that there remain linear combinations of 
F{km, kn) with numerical coefficients and rational additions. Indeed, by writing 

^F{km, kn) = °"' ~°" F(fe^, kn) + F{km, kn) (A.6) 
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and commuting (D^ — On) with the derivatives d/d\^ in (A. 2), one obtains 

°"^~°" F(0,fcO = (-l)^" ^^"~^^' +fenF(0,fcn-l) , A:„7^0 , (A.7) 
F{km,kn) — —kn — F[km — 1, kn) -\- {k^ + kn)F[km, kn — 1) (A. 8) 

+ kraikm + ~ ^)F{km ~ 1; ~ 1) ) 7^ , /c^ 7^ 

The latter identity apphed fc^ times with a subsequent use of (A.7) and (A. 3) brings to 
the following reduction formula: 

F{kmi kn) — {km'\~kn)F{kra-, kn + ( 1) " jz; , km 7^ , kn ^ 

(A.9) 



This formula is useful also in the case where there is a dimensional coefficient in 
front of F{krmkn) [18] since it enables one to convert the l/Dn into the l/Dm and vice 
versa. For that, it suffices to multiply Eqs. (A.7) and (A.9) by l/n^ ■ 

The relations above make it possible to bring the asymptotic form factors in [18] to 
their final form in (2.39)-(2.45). The strategy of this calculation is as follows. First, the 
relations (A.6)-(A.9) are used to get rid of the coefficients n^/n^ in front of F{km,kn) 
and, if needed, to replace the coefficients l/^n with l/D^ . Next, the relation (A. 3) 
is used in the thus obtained linear combinations of F{km, kn) to reduce the difference 
between km and kn ■ For example, 

F(0,4) + F(4,0) = 2F(2,2)-48F(1,1) + 24F(0,0) . (A.IO) 

The equality of the exponents km and kn ensures ultimately the positive definiteness of 
the total radiation energy. 

For illustration, consider two specific examples of the functions B^'' in [18]: 

and 

1 1 1 1 19 

B7\a2, ns) = ^F{1, 0) + -F(0, 1) + -F(2, 0) + -F(0, 2) + -F(l, 1) (A.12) 

1 1 n-^ / 1 1 X Do / 1 1 X 

F(3, 1) F(l, 3) + — (—F(2, 0) F(l, 1)) + — ( F(l, 1) + —F(0, 0)) . 

24 ^ ' ^ 24 ^ ' ^ ^2 ^96 ^ ' ^ 96 ^ ' Dg V 96 ^ ' ^ 96 ^ ' 
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In (A. 11), make the factors l/D in front of the F 's uniform, and use (A. 3). As a result, 
Bi2^ proves to be a tree operator: 

In (A. 12), first remove the factors □3/n2 , ^2/^3 with the aid of (A. 7) and (A. 9), and 
next use (A. 3). The result is 

S7^(n2,n3) = ^F(2,2)-lF(l,l) . (A.14) 

Apart from the simplification of the form factors, the expression for the effective action 
in [18] needs to be improved in the terms with the commutator curvature. In the work [IS- 
IS], the basis of nonlocal gravitational invariants was expressed from the outset through 
the source of the metric curvature by the use of the Bianchi identities. It is an omission 
of this work that a similar procedure has not been applied to the commutator curvature. 
By differentiating and contracting the Jacobi identities (1.5) one obtains the equation 

= y^ypn^-, - y^yanf,^ (a.is) 

which can be solved iteratively to express the commutator curvature TZajj through its 
source — ^'^T^a-y and the initial data. The iteration procedure is started by commuting 
the covariant derivatives on the right-hand side of (A. 15) with the aid of Eq. (1.10). 
Neglecting the commutators, one obtains to lowest order in 3? 

= V^jJc- V„J;3 + 0[5R2] . (A.16) 

The solution of this equation with the retarded Green function: 

7ta/3 = ^(V/jia - ^^Jp) + 0[3?'] (A.17) 

is the solution with no incoming wave of the vector connection field. (See Ref.[23] for the 
proof of a similar assertion in the case of the gravitational field.) 

From the results in [18] (with the vertex operators treated as above), the contribution 
of the commutator curvature to the scalar Eq. (2.39), obtains originally in the form 

^/(□i, D2){(^ + ^)v„7^l"^v^7^2^ + 7^l"^7^2«^} (a.is) 
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where 



f{a,,n,)^^F{2,2)-^F{l,l) . (A.19) 



The form (A. 18) is redundant because the commutator curvature is constrained by the 
Jacobi identities. By using (A. 17), one finds 

nr^T^2 a/3 = 2^ V^Ji°V^J2 a - 2^ J^^V^ j/ (A.20) 

where in the first term use can be made of the relation [15] 

2ViV2 = (Vi + V2)^-ni-n2 . (A.21) 

The contribution of the operator (Vi + V2)^ to (A.20) is of the form 

^(V, + V,fJ,°A,^n[(i>)(ii,)] . (A.22) 

Since here appears an overall □ operator acting at the observation point, the respective 
contribution to (A. 18) can be discarded (see Sec.2). As a result, expression (A. 18) takes 
the form 

-V;3Ji"V«J2'' . (A.23) 



□in2 

Finally, after the scalar I2 in (2.39) has been expressed as a quadratic combination of 
the conserved currents, the derivatives acting on the individual currents can all be made 
overall owing to the conservation laws (1.9). Thus, up to 0[K^] 

{W/3J,^){WJ,^) = WpWa{j{'J2^) , (A.24) 

(V«V;3 Jr) (V^V. J2"'') = VaV/5V^V.(jr^2"0 • (A.25) 
In this way the final result (2.39)-(2.45) is obtained. 
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Appendix B. Alternative approach to the non-scalar 

vertices. 



There is an alternative to the special consideration of the non-scalar vertices in Secs.4 
and 7. It is based on the fact that up to terms 0[^] the overall derivatives in (2.39) can 
be commuted with the operator log(— □) . Retaining only the vertex terms we have from 
(2.36) and (2.39) 

1 



rplXV 

vac 



(4^)2tr|V'^V'^l0g(-n)y.ealar 



+ V'^V'^V.V^log(-n)(cfss + Kl) 

+ V'^V''V„V^V^V,log(-n)T/^fr + trees| . (B.l) 

The derivatives act now at a point at X+ and, moreover, we need only the senior terms 
of expansions like (4.59) and (4.76): 

V^V^ ^ V^uV'u^ , V^V^V^Vfl ^ V^'uVuVauVau^ , etc. (B.2) 

Thus, making use of Eqs. (1.36) and (4.18), we obtain 

M(-oo)-M(oo) = -Jim^ /c?'5(0)tr ||^Dl(M,0|V;ealar) 

+ —Dr(u,(l>\Vfj'')VauVpuV^uV,u\ . (B.3) 

The price for getting rid of the derivatives is that now we have to deal with the moments 
Di of the non-scalar vertices. Since these vertices have different powers of growth at i"*", 
Eq. (6.41) is to be apphed in each case separately. Let us introduce a notation for the 
tensors at i+ that figure in the asymptotic expressions for the vertices, Eqs. (7.40)- (7.42): 

— OO 

— OO 
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180 



I / dT(^D-^)(^D^^) . (B.6) 



With the behaviours of the V 's in (7.40)-(7.42), the algorithm (6.41) yields 

4 1 



X, 



X I^[u, 0] 

X i+[7,0] 



for both V;?oi and , and 



n. — >no TzTT ^ J 



X (^"^(x, x)g%ix, x) A^^^^l 



a; 



X ^ [m, 0] 
X ^ i+[7,0] 



for V^grf^*^ ■ Here the appearance of 



(B.7) 



(B.8) 



(B.9) 



is owing to Eq. (4.16), and we didn't forget the propagators of parallel transport connQciing 
the point at X+ with the point at i'^ . 

The powers of u in expressions (B.7) and (B.8) are just the ones needed for the quantity 
(B.3) to be finite and nonvanishing but the weight of calculation transfers now to carrying 
out the integration over the directions of radiation i.e. over the angles at X+ . As seen 
from (B.3) and (B.7)-(B.8), one is to do the integrals 



d^S{(i)) (1 - 7nn)-5 
X (Vau{x)V(3u{x)g%{x,x)g'^^ 



X, X 



= [d'S{<f>){l-^nn)-' 

i+ J 



X — >• X+[m, 0] 
X i+[7,0] 



(B.IO) 



X 



(Vau{x)Vpu{x)Vju{x)Vau{x)g\{x, x)g^^{x, x)g'^^{x, x)g''^{x, x)) 



X 



X 



(B.ll) 

i+[7,0] 
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Since both points of the propagators g^ix, x) are at the future infinity, one may use the 
flat-spacetime expressions for these propagators. The integration is conveniently carried 
out in the Minkowski frame (1-25) for both points x and x . In this frame, 



g\{x,x)^S\ , Ve.«(x) = (5° - 



rii 



(B.12) 



and there emerge only integrals of the form 



J (f'S{(j)) (1 — 7nn) ^rii ■ ■ - rij 



(B.13) 



summarized in the table below. 

When contracting the H 's obtained in the Minkowski frame with the A 's in Eqs. 
(B.4)-(B.6), it should be taken into account that the Minkowski components of any of the 
A 's can be expressed as follows: 



A°- = 7V«r A^ 



(B.14) 



by virtue of the conservation equations (7. 23), (7.24). 
In this way after many cancellations we obtain 

47r 1 



IlapA-^ix) 

for both A^^'^gg and A"^t , and 
U^f,,,A'^^^'^{x 



3 (1-7^)3 



gapA-'^ix) 



:b.i5) 



1 



rg^f,g^,A^''^^'^\x) 



(B.16) 



5 (1-72)4' 

for A'^^J^^ where the latter tensor appears symmetrized. The use of these contractions in 
(B.7),(B.8) and (B.3) yields the same final result as in the main text, Eq. (8.4). 

The integrals (B.13) are elementary but it is useful to have them. Therefore, we 
present the ones needed for the calculation above. Denoting 



wy = ni(0)7* , 7i = Sik^'' , 7^ = 5ife7*7* 



we have 



(1 — 717)' 



(1 - 7^)^ 
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47r 
1 

An 



(fS{(j)) 

(fS{<i)) 



1 

47r 



1 

47r 



(1 


-717)4 




1 


(1 


— nj)^ 




1 


(1 


— nj)^ 




1 


(1 


— n^y 






(1 


— W))^ 




UiUk 


(1 


— n-)Y 




rii 


(1 


— n-^y 




TliUk 


(1 


— n-^y 


riiUkUp 


(1 


— ri'yy 


ninjUkUp 


(1 


— ri'jy 



1 (3 + 7^) 
3(1-72)3 ' 

(1 + 7^) 
(1-7^)4 ' 
1(5 + 1072 + 74) 

5' 



(1-^2)5 

1 (l + 372)(3 4 
3 (1-7^)6 

1 (5 + 7^) 
3(1-72)4^^ ' 

1 1 /r 7i7fc \ 

1 (35 + 4272 + 37^^ 

15 (1-72)6 
1 (5 + 37^) 8 (7 + 37^ 



-7i 



+ 



15 (1 - 


- 7^)'^ 


8 (9 


+ 7^) 


15(1 - 


_ y)6 


16 


1 


3 (1- 




8 


1 



7i7fe7p + 



15(1-72)6 

1 (7 + 7') 



15 (1 - 72)5 



15(1-72)5 

liljlklp + -TTTs r^^^ij^kp + ^ikSjp + 5ip5jk) 



15(1-72 

{kjlklp + ^ik^jlp + Klilk + ^jklilp + (^jp7i7fc + hplilj)- 



118 



References 



1. B.S.DeWitt, Phys. Rep. 19 C (1975) 295. 

2. A.A.Grib, S.G.Mamayev, and V.M.Mostepanenko, Quantum effects in intense ex- 
ternal fields (Atomizdat, Moscow, 1980), in Russian. 

3. N.D.Birrell and P.C.W.Davies, Quantum fields in curved space (Cambridge U.P., 
Cambridge, 1982). 

4. A.G.Mirzabekian and G.A.Vilkovisky, Phys. Lett. B 414 (1997) 123. 

5. V.P.Frolov and G.A.Vilkovisky, Phys. Lett. B 106 (1981) 307; in: Proc. 2nd 
Seminar on Quantum Gravity (Moscow, 1981), eds. M.A.Markov and P.C.West 
(Plenum, London, 1983) p. 267. 

6. J.Schwinger, in: Proc. 1960 Brandeis Summer School (Brandeis U., Brandeis, 1960) 
p. 282; J. Math. Phys. 2 (1961) 407. 

7. L.V.Keldysh, Zh. Eksp. Teor. Fiz. 47 (1964) 1515. 

8. Yu.A.Golfand, Yad. Fiz. 8 (1968) 600. 

9. P.Hajicek, in: Proc. 2nd Marcel Grossmann Meeting on General Relativity (Trieste, 
1979), ed. R.Ruffini (North Holland, Amsterdam, 1982) p. 483. 

10. E.S.Fradkin and D.M.Gitman, Fortschr. der Phys. 29 (1981) 381. 

11. J.L.Buchbinder, E.S.Fradkin, and D.M.Gitman, Fortschr. der Phys. 29 (1981) 187. 

12. R.D.Jordan, Phys. Rev. D 33 (1986) 44. 

13. E.Calzetta and B.L.Hu, Phys. Rev. D 35 (1987) 495. 

14. A.O.Barvinsky and G.A.Vilkovisky, Nucl. Phys. B 282 (1987) 163. 



119 



15. A.O.Barvinsky and G.A.Vilkovisky, Nucl. Phys. B 333 (1990) 471; Nucl. Phys. B 
333 (1990) 512. 

16. G.A.Vilkovisky, CERN-TH.6392/92; Publ. Inst. Rech. Math. Avancee, R.C.P. 25, 
Vol. 43 (Strasbourg, 1992) p. 203. 

17. A.O.Barvinsky, Yu.V.Gusev, G.A.Vilkovisky, and V.V.Zhytnikov, Covariant per- 
turbation theory (IV). Third order in the curvature (U. Manitoba, Winnipeg, 1993) 
pp. 1-192; J. Math. Phys. 35 (1994) 3543; Nucl. Phys. B 439 (1995) 561. 

18. A.O.Barvinsky, Yu.V.Gusev, V.V.Zhytnikov, and G.A.Vilkovisky, Class. Quantum 
Grav. 12 (1995) 2157. 

19. G.A.Vilkovisky, Class. Quantum Grav. 9 (1992) 895. 

20. A.O.Barvinsky, Yu.V.Gusev, G.A.Vilkovisky, and V.V.Zhytnikov, J. Math. Phys. 
35 (1994) 3525. 

21. A.G.Mirzabekian and G.A.Vilkovisky, Phys. Lett. B 317 (1993) 517. 

22. A.G.Mirzabekian, Zh. Eksp. Teor. Fiz. 106 (1994) 5 [Engl, trans.: JETP 79 (1994) 
1]- 

23. A.G.Mirzabekian and G.A.Vilkovisky, Phys. Rev. Lett. 75 (1995) 3974; Class. 
Quantum Grav. 12 (1995) 2173. 

24. A.O.Barvinsky and G.A.Vilkovisky, Phys. Rep. 119 (1985) 1. 

25. R.M.Wald, General relativity (Chicago U.P., Chicago, 1984). 

26. H.Bondi, M.G.J, van der Burg, and A.W.K.Metzner, Proc. R. Soc. A 269 (1962) 
21. 

27. R.Sachs, in: Relativity, Groups and Topology (1963 Lcs Houches Lectures), eds. 
C.DeWitt and B.S.DeWitt (Gordon and Breach, New York, 1964). 

28. A.M.Polyakov, Phys. Lett. B 103 (1981) 207. 

120 



29. J.L.Synge, Relativity: the general theory (North Holland, Amsterdam, 1960). 

30. B.S.DeWitt, Dynamical theory of groups and fields (Gordon and Breach, New York, 
1965). 

31. A.G.Mirzabekian, G.A.Vilkovisky, and V.V.Zhytnikov, Phys. Lett. B 369 (1996) 
215. 

32. E.S.Pradkin and G.A.Vilkovisky, Phys. Lett. B 73 (1978) 209. 

33. I.E.Tamm, Foundations of electricity theory (Gostekhizdat, Moscow, 1956), in Rus- 
sian. 

34. G.A.Vilkovisky, in: Quantum Theory of Gravity, ed. S.M.Christensen (Hilger, Bris- 
tol, 1984) p. 169. 

35. A.A.Ostrovsky and G.A.Vilkovisky, J. Math. Phys. 29 (1988) 702. 

36. N.M.J.Woodhouse, Phys. Rev. Lett. 36 (1976) 999. 



121 



Captions to the figures 



Fig.l. The Bondi-Sachs frame with the central geodesic r — and the hght cones of equal 
retarded time. 

Fig.2. Shown are the past light cone of a point x, an arbitrary point o on its surface, and 
both sheets of the light cone of o. (The line r = is the central geodesic.) The past 
light cone of x lies outside the light cone of o. 

Fig.3. The past hyperboloids of x, a{x, x) = q, inscribed in the past light cone of x, 

a{x, x) = 0. 

Fig.4. The past light cone of x in Fig.2 after x has moved to X+ along the radial geodesic 
shown bold. The resultant null hyperplane has the parameters u, 4> of this geodesic. 

• 

Fig.5. Shown are the support tube of a physical source J, the compact support of J, and 
the earliest hyperboloid of x that still crosses the latter support. The broken lines 
bound the causal future of the support of J. 

Fig. 6. Penrose diagram for the section of fixed angles in the Bondi-Sachs coordinates. The 
exterior of the light cone of o is divided into four subdomains bounded by two future 
light cones u — Ui and u — U2, and the surface of the tube r < Tq. 

Fig. 7. Lorentzian section F of a spherically symmetric spacctimc. The broken line is the 
mapping on F of a timelike radial geodesic. The lines T"^ — r bound the mapping 
on F of a spacehke hyperplane. The bold hne is the world hue of the support shell 
(or the boundary of the support tube) of the source. 
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Fig.2 
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Fig.4 
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Fig.5 
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Fig.6 
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